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The critical behavior of semi-infinite d-dimensional systems with n-component order parameter 
(j} and short-range interactions is investigated at an m-axial bulk Lifshitz point whose wave-vector 
instability is isotropic in an m-dimensional subspace of R''. The associated m modulation axes are 
presumed to be parallel to the surface, where < m < d — 1. An appropriate semi-infinite |</)|^ 
model representing the corresponding universality classes of surface critical behavior is introduced. 
It is shown that the usual 0{n) symmetric boundary term oc </)^ of the Hamiltonian must be sup- 
plemented by one of the form A ^^^i{d4> / dxa)^ involving a dimensionless (renormalized) coupling 
constant A. The implied boundary conditions are given, and the general form of the field-theoretic 
renormalization of the model below the upper critical dimension d*{jn) — 4-f m/2 is clarified. Fixed 
points describing the ordinary, special, and extraordinary transitions are identified and shown to be 
located at a nontrivial value A* if e = d* (m) — d > 0. The surface critical exponents of the ordinary 
transition are determined to second order in e. Extrapolations of these e expansions yield values 
of these exponents for d = 3 in good agreement with recent Monte Carlo results for the case of a 
uniaxial (m = 1) Lifshitz point. The scaling dimension of the surface energy density is shown to be 
given exactly hy d + m{0 — 1), where 9 = 1^14/^12 is the anisotropy exponent. 
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I. INTRODUCTION 

As is well known, there are many physical systems 
whose phase diagrams exhibit a Lifshitz point. A Lif- 
shitz point is a multicritical point at which a disordered, 
a homogeneous ordered, and a modulated ordered phase 
meet. It divides the critical line that separates the disor- 
dered phase from the two ordered ones into two sections, 
and marks the onset of a wave-vector instability. In the 
case of an m-axial Lifshitz point this instability occurs in 
an m-dimensional subspace of M'*. Prominent examples 
of systems with uniaxial (m = 1) Lifshitz points are, on 
the experimental side, the metallic compound^"^ MnP, 
and on the theoretical side, the three-dimensional axial 
next-nearest-neighbor Ising (ANNNI) model. ^'^^ 

Although the notion of a Lifshitz point was introduced 
decades ago,^^ the successful application of modern field 
theory tools to the study of critical behavior at such 
points is a fairly recent development. ^^"^^ A full two- 
loop renormalization group analysis and e expansion of 
all critical exponents to second order in e about the up- 
per critical dimension d*{m) = 4: + m/2 (0 < m < 8) has 
been accomplished for general values of m only two years 



ago 



14,15 



In the present paper we wish to study the effects 
of surfaces on the critical behavior at an m-axial Lif- 
shitz point. Previous work on this problem is scarce. 
It started in 1986 with Gumbs' investigation^^ based 
on Landau theory. This was reconsidered, completed, 
and partly corrected by Binder and Frisch.-'^^ In either 
one of these papers, only the uniaxial case m = 1 
was considered, and the special axis along which ferro- 
magnetic nearest-neighbor and antiferromagnetic next- 



nearest-neighbor interactions compete was chosen per- 
pendicular to the surface plane. The complementary case 
in which the special axis of the considered semi-infinite 
ANNNI model is parallel to the surface was dealt with 
on the level of mean field theory in subsequent work 
by Frisch et al?^ Even more recently, Pleimling^^ pre- 
sented results of Monte Carlo simulations of semi-infinite 
ANNNI models for both types of surface orientations. 

A characteristic feature of critical behavior at Lifshitz 
points is anisotropic scale invariance. Let Axa and Axp 
with 1 < a < m and m < (3 < d denote displacements 
along the corresponding Euclidean axes and suppose that 
the Aa:^ are rescaled by a factor £. Anisotropic scale 
invariance means that the Axa must be rescaled by a 
nontrivial power £^ in order that the system looks self- 
similar on different scales. 

An obvious consequence of this property is that two 
distinct cases of surface orientations must be distin- 
guished when dealing with surface critical behavior at 
bulk Lifshitz points: a parallel one for which the surface 
normal n is along a /3-direction (all a-directions being 
parallel to the surface plane), and a perpendicular one 
for which n is along an a-direction. Since the distance z 
from the surface scales differently in these two cases, one 
expects (i) that the respective values of the surface expo- 
nents are in general different, even in mean-field theory, 
and (ii) that the boundary terms which must be included 
in the Hamiltonian are different as well. Expectation (i) 
is borne out by both the mean-field results of Refs. 19 
and 20 as well as by Pleimling's Monte Carlo results;^^ 
(ii) turns out to be equally true.^'^ 

In this paper^^'^^ we will focus our attention on the 
case of parallel surface orientation, leaving the study of 
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the other, "perpendicular" case to a subsequent paper 
[which will confirm the above conclusion (ii)]. To our 
knowledge, experimental investigations of even bulk crit- 
ical behavior at Lifshitz points^'^ have so far been re- 
stricted to the uniaxial case m = 1, and this is also still 
the only one for which Monte Carlo simulation results 
on surface critical behavior at bulk Lifshitz points are 
available. 

Although m = 1 thus appears to be the case of great- 
est interest, we will keep the value of m as general as 
possible, imposing at this stage no restriction other than 
the obvious one < m < d ~ 1.^^ Just as in essentially 
all renormalization group (RG) analyses of bulk critical 
behavior with m > 1, we will make, however, the sim- 
plifying assumption that the wave-vector instability is 
isotropic in the m-dimensional subspace of a directions. 
Let us nevertheless add a cautionary remark about this 
presumed "m-isotropy" . In a recent study of bulk criti- 
cal behavior, this assumption has been relaxed by re- 
quiring only invariance in this subspace under the cubic, 
or even a lower, symmetry group. The results indicate 
that the isotropic bulk fixed point for m > 2 gets desta- 
bilized by bulk terms of second order in <p and fourth 
order in the derivatives d/dxa, although the associated 
crossover exponent appears to be fairly small. The in- 
clusion of such non-isotropic terms is beyond the scope 
of the present work and, of course, meaningless in the 
uniaxial case m = 1. 

The remainder of this paper is organized as follows. In 

the next section wc specify the Hamiltonian representing 
the corresponding universality classes for surface criti- 
cal behavior at m-axial bulk Lifshitz points. We explain 
which boundary terms must be included in it and give the 
implied boundary conditions. In Sec. Ill we first clarify 
the renormalization of the model for general values of its 
parameters and derive the resulting RG equations. On 
the basis of two-loop results, we then discuss the form 
of the RG flow in the space of the surface interaction 
constants and identify the fixed points describing the or- 
dinary, special, and extraordinary transitions. 

Section IV deals specifically with the ordinary transi- 
tion. We show that the RG analysis can be simplified 
in much the same way as in the m = case of a criti- 
cal point^^ by choosing Dirichlet bomidary conditions, 
considering correlation functions involving the normal 
derivative of the order parameter at the boundary, and 
making use of the boundary operator expansion. The re- 
quired RG functions are evaluated to two-loop order; the 
scaling index of the surface energy density is determined 
exactly. 

In Sec. V results to order are given for the critical 
exponents of the ordinary transition. These e expansions 
are exploited to estimate the values of these exponents 
for the uniaxial, one-component case m = n = 1 in three 
dimensions. Section VI contains a brief summary and 
concluding remarks. Finally, there are five appendixes, 
explaining details of our calculations. 



II. THE HAMILTONIAN, ITS BOUNDARY 
TERMS, AND IMPLIED BOUNDARY 
CONDITIONS 

We consider systems with short-range interactions, as- 
suming also that the perturbations of these interactions 
induced by the presence of the surface decay to zero 
within a short distance from it. With these assumptions 
the Hamiltonian can be taken to be of the form 

H = / £b{x)dV+ [ £i{x)dA, (1) 

where Ch{x) and Ci{x) depend on the order parame- 
ter density 4>{x) = {(j)a{x), a = l,...,n) and its spa- 
tial derivatives up to a finite order. Here the volume 
and surface integrals extend over 2J, the d-dimcnsional 
half-space = IR'^-ix[0,oo), and 03, the d - 1 di- 
mensional surface plane z = 0, respectively. Writing 
X = ((.Xq), (a;^)) = (fTz), we split the position vec- 
tor X into its m-dimensional component (xa) € 
and {d — m)-dimensional one (x^) = ((r/3),2;), where 
= ((^a) = ixa),{rp)) is the (d — l)-dimensional co- 
ordinate along the surface. We choose the same bulk 
density as in Refs. 14-17, namely 

o / m \ 2 -. d 

o m o o 

a=l 

where da and dp denote the spatial derivatives d/dxa and 
d/dxf) with 1 <a <m and m+1 < < d, respectively. 
Had we not assumed that the wave-vector instability is 

isotropic in the subspace M™, the term {Yl^=\^'a4') 
would have to be supplemented by similar, albeit less 
symmetric, terms involving four derivatives, such as^'^ 

In order to decide which monomials should be included 

in Ci{x), we use power counting. Recalling from Ref. 14 
the naive dimensions [xfj] = [2] = [xa\ = ct^^'* A*"""^^^, 
and [(j){x)] = (T-^/s ^(d-2-m/2)/2 (^here /x is an arbi- 
trary momentum scale), we include only such 0{n) in- 
variant monomials whose interaction constants have non- 
negative /i-dimensions at the upper critical dimension 
d*{m). It is not diflacult to see that the choice 

A(a.) = |0^ + |Ea0f (3) 

a=l 

is sufficient since the two additional candidates ^dada<j> 

and (pdnCp both may be dropped. The integral of the 
former reduces to that of the derivative term retained in 
Eq. (3) upon integration by parts, and the latter can be 
shown to be redundant in much the same way as in the fa- 
miliar case of the standard semi-infinite 4>* model. One 
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must merely note that the action defined by Eqs. (l)-(3) 
leads to the boundary condition^*"'^^ 



dn4>{x) = (c — A dada) 4>{x) , X £ ^ 



(4) 



which in turn implies that the surface term (pdncj) 
is equivalent to a linear combination of the two mono- 
mials present in Ci. Here a sum convention which 
will frequently be employed below has been introduced: 
Pairs of equal a and (3 indices are to be summed over 
a = 1, . . . ,m and /? = m + 1, . . . , d, respectively. 

We need the free propagator of the disordered phase, 
G{x,x') = [5'^n/d<p6<p\^=o\~Hx^x'). This is a solution 
to the equation 

[a {dadaf - P dada ~ dpdfj + f] G{x, x') = 6{x - x') , 

(5) 

subject to the boundary condition 

{dn-C+°Xdada)G{x,x')=0 , X € ^ , x' ^ ^ , 

(6) 

where dn (= dz) denotes the derivative along the inward 
normal n. 

We take periodic boundary conditions along all — 1 
axes parallel to the surface. Let p = {{pa), (pp)) G 
j^m^-j^d-rn-i ^j^g momentum conjugate to r, and 
G{p; z, z') be the corresponding Fourier transform of G 
with respect to the c? — 1 coordinates parallel to the sur- 
face. In this pz representation G is easily evaluated; one 
finds 



G{p; z, z' 
1 



-Kp|2 



-z'l _ £p_ 



with 



and 



Cnp 



C+\paPo 



^JP g-Kp(z-|-Z ) 



PPaPa +P[jP[j 



O (paPa 



(7) 



(8) 



(9) 



The part of G that depends on \z — z'\ is its bulk ana- 
log Gb; the remainder is the contribution induced by the 
presence of the surface. Results for the free bulk propa- 
gator Gh{x — x') in position and momentum space can be 
found in Refs. 14 and 15. Specifically at the (Gaussian) 
Lifshitz point f = p = 0, one has 

with 



(10) 

(11) 



and 



^m,d(w) 22-l-m7r(6-|-m-2e)/4 
( 



X 1F2 



^ e 1 m + 2 v'^ 
~ 2' 2' 4 '64 



.'5-f 
2 



2 '2' 



T' 64/ 



4r(i + f) 

(12) 



where e = 4 + ^ — d while pFg means the generalized 
hypergeometric function. 

In our calculations below we need, in particular, the 
free bulk propagator between a point x = (r, z) and its 
mirror point x = (r, —z) = x — 2zn. With the aid of the 
Taylor expansion of the scaling function ^m,d{y) given in 
Eqs. (10) and (11) of Ref. 15 one easily obtains the result 



Gb{x -x) = Gb{2zn) = $„,d(0) a'^/^ (20)^-^ 
^ r(2-e) sin(e7r/2) , 2 



(13) 



in which 



p _ r(i + 6/2)r^(i-./2)r(m/4) 

(4^)(8+m-2.)/4r(2-e)r(m/2) ^ ' 

is a factor introduced in Ref. 15; just as there, \vc will 
absorb it in the renormalized coupling constant it to be 
defined below. 



III. RENORMALIZATION GROUP 

A. General considerations 

Let us introduce the CTimulants involving N fields 
3^ at points Xj off the surface and M boundary fields 



^H^'^k) = 06fc(rfe,O), namely 



N 



M 



GW^)(x;r) = (n<^a,(^i)n€(^fc) 



fe=i 



(15) 



Here x and r are convenient short-hands for the sets of 
position variables {xj} and {rk}, and the components 
indices {aj} and {bk} have been suppressed on the right- 
hand side for the sake of notational simplicity. 

We wish to analyze the critical behavior of these func- 
tions in bulk dimensions d < d* (m) using a field-theoretic 
RG approach and the e = d* (m) — d expansion. To reg- 
ularize their ultraviolet (uv) singularities, we employ di- 
mensional regularization. Aside from the bulk uv singu- 
larities induced by the bulk part of the free propagator, 
additional primitive ones localized on the surface occur. 

Consider first the former "bulk uv singularities" . In 
previous investigations^^' ^'"^ of the bulk model, p and 
T were set to their values plp and tlp at the Lif- 
shitz point in the actual calculations, or else deviations 
df = f — Tlp ^ with 5p = p — p^p = and 5f ^ with 
6p = were considered. This is sufficient to determine 
the two RG eigenexponcnts 1/j^;2 and (f/i'12 associated 
with the corresponding cigcnopcrators at the infrared (ir) 
stable fixed point (where 1^12 is a standard correlation ex- 
ponent while (p means the crossover exponent pertaining 
to 6p). One must remember, however, that in the p- 
dependent scaling forms of the correlation functions de- 
rived in Refs. 14 and 15 the linear scaling fields associated 
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with these cigcncxponents ought to be replaced by non- 
hnear ones when considering general deviations S 'p and 
5f from the Lifshitz point. 

Here we wish to go beyond these previous analyses 
by allowing both and Sp to be nonzero and not nec- 
essarily small. Since the naive dimension of pa~^l'^ is 
\pb~^l'^\ = fj,, we see that deviations Sp may give 
contributions to the renormalization of ^f. That is, 
counterterms of the form oc /ip^ 0^ (where p is the 
dimensionless renormalized counterpart of 5pa~^/'^) are 
possible. Combining this with the considerations made 
in Ref. 14 and 15, we can conclude that the bulk uv sin- 
gularities can absorbed by making the following "bulk 
reparametrizations" : 





yl/2 , 

~ ^4> V'ren ) 




= Z^a , 


T - TLP 


= 1? Zr [t + Ar p- 


(p-pLP)a-i/' 


= P-ZpP , 







(16) 



Here Z^ = Z^{u, e), l = a, p, u are bulk renormalization 
factors for which results to two-loop order were given in 
Refs. 14 and 15 for general values of m. In our pertur- 
bative approach based on dimensional regularization and 
the e expansion the critical values tlp and pLP vanish, 
as usual. At is the additional renormalization function 
associated with the counterterm discussed above. Owing 
to the restrictive assumptions made in Ref. 14 and 15, it 
was not needed — and hence not computed — there. A 
straightforward one-loop calculation yields 



Ar{u,e) = ^ +0^2 

3 16 e 



(17) 



Next, we turn to the "surface uv singularities" in- 
duced by the free propagator's "surface part", i.e., the 
part proportional to e^'^p'^^^ L They require addi- 
tional counterterms, localized on the surface. To deter- 
mine their form, we use power counting in conjunction 
with what is generally known about such field theories 
with planar boundaries. ^^'^^ In comparison to the usual 
(m = 0) semi-infinite (/)^ model, a qualitative change oc- 
curs: There is a surface variable with zero /^-dimension, 
namely Air^^/^. This implies that the surface renormal- 
ization functions do not only depend on the renormalized 
bulk coupling constant u but additionally on the renor- 
malized analog of this bare variable, which we denote as 
A. Furthermore, power counting suggests that a surface 
counterterm of the form oc dA J2a('^<^^)^ needed 
even if A is set to zero. If a large-momentum cutoff A 
were used to regularize the uv singularities of the the- 
ory, then the associated renormalization function would 
diverge as In A. Hence it must have poles in e in renor- 
malization schemes based on dimensional regularization 
such as ours. 



Finally, arguments completely analogous to those giv- 
ing the contributions (x to the renormalization of f tell 
us that the renormalization of the surface enhancement 
variable c involves contributions linear in p, and hence is 
not multiplicative. The upshot of these considerations is 
that the following reparametrizations of surface quanti- 
ties are needed: 



A<7-i/2 



p,Zc [c-\- Ac{u,X,e) p] 
\ + Px{u,X,e) . 



(18) 



Here the surface renormalization factors 
Zi^c{u, A, e) have the form 



Zi, 



We 



oo oo 
i,j=l k=0 



provided we fix them by requiring that the poles be min- 
imally subtracted. Likewise we have for the renormaliza- 
tion function Px, 



oo oo 



i,j=l k=0 



and similarly for Ac- 

It is easy to see that the one-loop coefficient pI^^'~^^ (A) 
vanishes for A = 0. This follows from the fact that the 
tadpole graph Q is independent of the momentum p, so 
that no subtraction corresponding to a surface two-point 
counterterm oc p'^ is required for the one-loop graph of 
{(t)(fy^) if A = 0. On the other hand, a subtrac tion of 
this kind is needed for its two-loop subgraph 
with A = (where the crossed circle marks the external 
surface point). Hence, ' (A = 0) does not vanish. 



B. One-loop results for general values of A 



At one- loop order, the bulk renormalization factors Z^, 
Zo-, Zt, Zp, and Z„ are known to be independent of m 
and hence equal to their m = analogs. Obviously, this 
property cannot be expected to carry over to the surface 
renormalization factors Zi {u, A, e) and Zc{u, A, e) because 
of their A dependence. On the other hand, it does hold 
for Zi(u,0,e) and Zc(u,0,e). 

To see this, note that a straightforward calculation pre- 
sented in Appendix C yields the one-loop results 



Zi{u,X,e) = 1 



n + 2 ii{X; m) u 
~3 2e 



+ 0(w2), (21) 



5 



ZM, A, e) = 1 + !i±^ [2^,{X■,m)-^.{\ra)\u ^ ^ 

(22) 

and 



P\{u, A,e) = - 



n + 2 H(A;m)A» ^ _ ^^3) 

3 2c 



Here ii(A; m) = ii(A, 0; m) and i2(A; m) = i2(A, 0; m) are 
special cases of the integrals 



/•! 2t("'-2)/2 (1 - /2) 

ii(A,e;TO) = dt 
Jo 



9> (2-2(-'„0/4 



1 - 

5[m/4, (6 - m - 2e)/4] 1 + At 

(24) 



and 



i2{X,e;m) 



1 2t(™-2)/2 (1 _t2^(2-2e-m)/4 ^ 

dt ■ 



S[m/4, (6 - m - 2e)/4] (1 + At) 



where B{o.. h) is the Euler Beta function. 
At A = and for A — > oo wc have 



2 ' 



(25) 



ii{0;m) = i2(0;m) = — ii(oo;m) = 1 , i2(oo;m) — . 

(26) 

These values for A = ensure that our one-loop results 
(21) and (22) for Zi and Z^ reduce to their m = 
analogs^^ when A is set to zero. 



C. RG equations 

Upon changing /i at fixed values of the bare inter- 
action constants, we see that the renormalized func- 

tions G^e^n^) = ^-(iV+M)/2^-M/2^(;V,M) ^^^isfy the 

RG equations 



N + M M 
+ — 2 — ^'^ + Y ^1 



GW^)=0 (27) 



with 



'Da = ndu. 



=u^a,r,p,c,X 



(28) 



Here r]^{u) and 771 (u, A) denote particular ones of the 
functions 



= <P,u,a,T,p, l,c, 



(29) 



where 1-18^1^ stands for a /j. derivative at fixed bare inter- 
action constants m, f, p, o", c, and A. 

The beta functions appearing in Eq. (28) are defined 
via 

/?p=/ia^|oP, p = U,(7,T,p,C,X. (30) 

They can be conveniently expressed in terms of the 77^, 
br{u) = Ar [ /xS^lglnA^ + i1t- 2r/p] , (31) 



and 



bc{u, X)=Ac[ p,d^i\Q In Ac + % - 'Hp] ■ (32) 



We have 



(3u{u,e) = -u[e + r]u{u)] , 

(3a{u,a) = -ariaiu) , 

l3r{u,T,p) = -T[2 + 'qr(.u)]- p^hr{u) , 

(3p{u,p) = -p[l + r]p{u)] , 

(3c{u,X,p,c) = -c[l + T]c{u,X)]- pbc{u,X) , (33) 



and 



0x{u,X) 



-f3u{u, e) duPxju, X, e) 
l + dxPx{u,X,e) 



(34) 



That the functions r}u,a,T,p{u), br{u), r7c,i('M, A), bc{u,X), 
and f3x{u,X) arc independent of e is due to our use of 
the minimal subtraction prescription. As usual, the eta 
functions can be written in terms of u derivatives of the 
residues of the Z factors; we have 

Vp{u, A) = -udu Res Zp{u, A, e) , p = u, u, r, p, 1, c . 

(35) 

The analogous results for b-r, be, and /3a, 
&x(w) = — u9„ Res y4.x(u, e) , 

e=0 



and 



bc(u,X) = — u9„ResAc(u, A,e) , 

e=0 



f3\{u,X) = uduResPx{u,X,e) 

e=0 



(36) 



(37) 



XzPf-^)(A), 



(38) 



can be derived from Eqs. (31), (32), and (34), respec- 
tively. 

Upon substituting Eq. (17) into Eq. (36), we obtain 



, , , n + 2 mu , ^, 2n 



(39) 



Our pcrturbative result for f3\ will be given and discussed 
below [sec Eqs. (55) and (59)-(61), and Fig. 1]. The 
function bc{u) will not be computed in this paper since 
its explicit form is not needed for our subsequent analysis. 



D. Flow equations and fixed points 

To exploit the RG equations (27) via characteristics, 
we introduce running coupling constants p{t) into which 
the p evolve under a change p p{i) — p£ of the mo- 
mentum scale. They are solutions to the flow equations 



i-^p{e) = f3p[u{£),...] , p = u,a,T,p,c,X, (40) 
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satisfying the initial conditions 

P(1) = P, P = U,(^,T,P,C,X. 



(41) 



As is well known, for dimensions d < d*{m) the bulk 
critical behavior at the Lifshitz point r = p = (with 
cr > 0) is described by a fixed point that is located at 
the nontrivial zero u* = u*{m, e) of the beta function 
According to Eq. (60) of Ref. 15, the expansion to order 

of this root is given by 



2e 



u = 



3 n + 8 
1 n + 
24 ~3 



8 £2 


9 


^ 5n + 22 


27 


_n + 8_ 


r 27 



+ 



3a{m) 
8 (m + 2) 



Ju{m) 



(42) 



where Ju{m) is one of the four single integrals j^{m), 
jai'iTT'), jpifTT-)! ^iiid Ju (m) in terms of which the two- loop 
results for the bulk renormalization factors Z^{u), Z„{u), 
Zp{u), Zt{u), and Z„(u) were expressed in Ref. 15. It is 
given by 



Ju{m) 



1 - 



+ ju{m), (43) 



where Ce = 0.577216 . . . denotes Eulcr's constant while 
ip{x) is the digamma function and jt,(m) means the in- 
tegral 



juim) = 
with 

and 



Br, 



poo 

Jo 

(44) 



32 r(i + 2?) ' H ' '2' '2^ 4 '647 



4 r(i + ^) 



^^H2'2'^+4'64)(^^^ 



Bm = 



TP2 



2W+m ^6+^ r(Y) 

r(2-f)r(f)2 



(46) 



The quantity 5^ denotes the surface area Sm = 
2 7r™/^/r(m/2) of an m-dimensional unit sphere. 

Since some of our analytical results to be given below 
involve besides the functions j<^(m) and jo- (w) alsojp(m), 
let us recall their definition here for completeness. We 
have 



and 



j^(m) = Jo,3(m) 



ja{m) = Bm ^4,3 (to) , 



(47) 



(48) 



where Jo,3(m) and J4,3(m) are special cases of the inte- 
gral 



Jp,s{m) 



oo 



(49) 



previously considered (even for d ^ d*) in Ref. 14. Fur- 
ther, 



Jp{m) = Bm / dvv"'+^ ^m,Av) ^mA'{v) , (50) 
Jo 



with 



iv) 



V 



(4-m)/2 



32 (27r)- 



/™-4 - 



I — 



(51) 

where Ii, and Ljy denote modified Bessel and Struve func- 
tions, respectively.''^ 

Here we are interested in the surface critical behavior 
at the analog of the ordinary transition that occurs at the 
Lifshitz point (LP). Whenever it is necessary to distin- 
guish this type of ordinary transition from its counterpart 
taking place on the ferromagnetic section of the critical 
line (CL), we shall refer to the former as LP ordinary 
transition and to the latter as CL ordinary transition. 
The latter belongs, of course, to the surface universality 
class of the usual ordinary transition at a critical point 
(CP), which is described by the standard semi- infinite 
n-vector model. 

Let us note some characteristic features of the LP ordi- 
nary transition, which serve to identify it and distinguish 
it from other types of surface transitions such as the spe- 
cial and extraordinary ones: 

(i) The surface is disordered on the bulk disordered 
side of the transition. Hence the order-parameter profile 
{(j){x)), and especially the local surface order parameter 
((/)®(r)) = {(p{r,z=0)), vanish at the transition. 

(ii) The surface does not become critical on its own 
but remains noncritical at the transition. That is, the 
ir singularities which occur at the surface are induced by 
the bulk, and the surface susceptibility 



(52) 



diverges neither at p = nor at any nonzero value of the 
momentum p, so that its inverse is strictly positive: 



Xri'(p)>0, Vp. 



(53) 



In the case of the LP special transition, condition (i) 
continues to apply. However, Eq. (53) holds only for 
nonzero momenta p because the surface becomes criti- 
cal as well; i.e., the surface susceptibility xii = Xii(O) 
diverges, XiiW = 0- O^i the other hand, at the LP 
extraordinary transition, mi = {(f>^) does not vanish. 
Therefore (i) ceases to hold while (ii) still applies for all 
p (barring eventual Goldstone singularities at p = due 
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-ii(X;m)X 



m 




FIG. 1: The functions — ii(A; m) A for m = 1, 2, . . . , 6. 



to the spontaneous breaking of the continuous 0{n) sym- 
metry for n > 2 in sufficiently high dimensions). 

From the form of the beta function /3c in Eq. (33) for 

p = we can read off the fixed-point values c* = ±oo 
and c* = of c; all others would require that 



77e(w*,A*) = -l 



(54) 



We ignore this possibility since we see no reason why 
the zeros A = A* of the function f3\{u* ,\) should satisfy 
Eq. (54). (By inspection of the perturbation series for rj^ 
and I3\ to low orders in u one can convince oneself that 
this condition cannot hold as an identity in e.) 

Recalling that the line of CP ordinary transitions is 
mapped onto a fixed point with c = d^^^ = oo, we expect 
the fixed point describing the critical behavior at the LP 
ordinary transition to be located at c = oo as well. In 
the special case m = 0, in which the variable A drops 
out, this is evident because the results must reduce to 
those for the usual isotropic semi-infinite cj)'^ model (cf. 
Refs. 28 and 36-39). However, for m > this remains to 
be verified. 

To this end we need some information about the flow 
in the cA hyperplane a.t u — u* and p = t = 0. Upon 
inserting the 0{u) result (23) for P\ into Eq. (38), we 
arrive at a two-loop expression of the form 



/Ja(u,A) = --^ii(A;m)Au + 2u2p(2>-i)(A)+0(w3) . 

(55) 

Let us first consider (3\ (u* , A) in the one- loop approxi- 
mation, i.e., to first order in u. Its fixed-point values are 
given by the zeros of the functions — Aii(A; m), which are 
plotted for m = 1, 2 . . . , 6 in Fig. 1. 

Prom their form we see that the trivial zeros A = are 

ir unstable in the A direction. There also exist nontrivial 
positive roots A!j_ for 1 < m < 6, given by the zeros Ao(to) 



of the functions ii(A; m): 



ii{Xo;m)=0 for Ao(rn) 



' 6.2092 


for 


m 


= 1, 


2.5129 


for 


m 


= 2, 


L7018 


for 


m 


= 3, 


1.3425 


for 


m. 


= 4 , 


1.1363 


for 


m 


= 5, 


. 1 


for 


m 


= 6 . 



(56) 

Since (3\{u* , X) has positive slopes at Ao(to), these fixed 
points are ir stable in the A direction. 

At two-loop order the trivial zeros get shifted to 



A* 



y2p(2,-l)(0) 

{n + 2){n + 8) 



(57) 



In deriving this result we have substituted u* by it e 
expansion (42). 

The associated RG eigenvalues (which govern the be- 
havior X(£) — A* ~ £~y^ of the running variable A near 

A,) are 



yx = -{9xPx){u*,X*) = 



n + 8 



e + 0{e^) . 



(58) 



According to this 0{e) result, > for e > 0. Thus the 
fixed points with A = A* are indeed ir unstable. 

The two-loop function P^^'~^\X) is as yet unknown, 
but in Appendix D we calculate its value at A = 0, show- 
ing that it is negative. In the special cases m = 2 and 
m = 6, it can be computed analytically. One obtains 



p(2,-l)(0) 



2 1 



TO=2 



3 192 

,2 



[57r2 - 16 In 2 



-61n^3- 12Li2(l/3)] (59) 



and 



P(2.-l)(0) 



n-l-2 1 
3~ 6 ' 



(60) 



where Li2(a;) = X^fcli^'^/^^ is the dilogarithm, giving 
Li2(l/3) = 0.366213 For other choices of m, the in- 
tegrals in terms of which P(^'~^^(0) is expressed in Ap- 
pendix D [see Eqs. (B5), (B6), and (D7)] can be deter- 
mined by numerical means (see Appendix E). Our results 



p(2,-l)(0) 



' 0.12473 


for 


m 


= 1, 


0.13865 


for 


m 


= 2, 


0.14885 


for 


m 


= 3, 


' 0.15652 


for 


m 


= 4, 


0.16231 


for 


m 


= 5, 


. 0.16667 for 


m 


= 6, 



(61) 



demonstrate that P(^'~^)(0) does not vanish. 

In order to understand the meaning of the fixed point 
at {u, c, A) = (m*, 0, A*), one should note that the deriva- 
tive dxii {p)/d{paPa)\p_Q of the zero-loop inverse sur- 
face susceptibility vanishes at A = 0. This behavior of 
Xii (p) is reminiscent of the vanishing of the derivative 
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df^^ {q, T=0, p) I d((iaqo)\ ^ ^ of its bulk analog, the bulk 
vertex function f[,^^(q) = \ d'^ xY^^^ (x) e^'i '^ , at p = 0. 



Clearly, if A becomes negative, then stabilizing surface 
contributions to the Hamiltonian of fourth order in da 
and second order in will be needed. In sufficiently high 
space dimensions d, a phase with modulated surface or- 
der and a surface Lifshitz point should exist for such a 
generalized model, at least when the bulk parameters t, 
/9, and u are in a regime for which the bulk is disordered. 
Beyond Landau theory, fluctuations arc expected to shift 
the location of the "instability point" of the surface sus- 
ceptibility (i.e., where 9xn"^(p)/9(PaPa)|p_Q vanishes if 
c = p = T = 0)toa nonzero value of A. Hence it is 
reasonable that the fixed point associated with this sur- 
face wave-vector instability gets shifted to the nontrivial 
value (57) of A. 

In d = 3 (bulk) dimensions, no surface Lifshitz point 
is expected to occur at absolute temperatures T > 
for models of the kind considered (with short-range in- 
teractions) because the bulk phase in d — 1 dimensions 
to which the required surface phase with long-range 
modulated order would correspond to is believed to be 
thermodynamically unstable. Thus the fixed point at 
(?i. A) = (m*, a*) and c = /9 = r = 0isof little interest. 
Since our main concern here is the ordinary LP transi- 
tion, we will refrain in the sequel from a detailed inves- 
tigation of the multicritical behavior it may describe in 
dimensions d > 3. 

Next, we consider the nontrivial zeros A^ of /3a(A, u*). 
Using the two-loop expression (55) in conjunction with 
the previously utilized result (42) for u*, one obtains 



A!;(m) = Ao-(-- 



72 p. 



(2,-1) 



A 



(Ao)e 



n-|-2)(n + 8)Ao «i(Ao;m) 



-O(e^) , (62) 




FIG. 2: Schematic picture of the RG flow in the cA plane at 
p = T = and u = M*, showing the fixed points Pprdi 



and specified in Eq. (65). 



73* 

' sp: 



All three arc ir stable along the A direction. As the nota- 
tion suggests, we therefore expect them to describe the 
LP ordinary, special, and extraordinary transitions, re- 
spectively. 



IV. RG ANALYSIS OF THE LP ORDINARY 
TRANSITION 

A.. Reduction to cumulants with c = 00 and A = 



where i\[\\rn) = dii{\;m) / d\. 

For the associated correction-to-scaling exponent 



(63) 



which governs the behavior A — A!j_ ~ of the running 
variable A near A!j_, we find the result 



i'i(Ao;m)Aoe-FO(e2) . 



(64) 



By performing the differentiation d/d\ inside of the 
integral (24), or from Fig. 1, one can easily see that 
i'i{X;m) < for A > 0. Hence wa > to linear order 
in e > 0. 

A schematic picture of the flow is depicted in Fig. 2. 
On the line A = in the hyperplane u = u* we can 
identify the three fixed points 



' ord 
sp * 



(c^p = 0, A = a;) , 
(c^x = -00, a = a;) . 



(65) 



Assuming that the initial values u, c, and A belong 
to the basin of attraction of the fixed point V*^^, we 
now specialize to the case of the LP ordinary transition. 
From the explicit form (7) of the free propagator G and 
the boundary condition (6) it satisfies we see that in the 
limit c ^ 00 at fixed A > 0, the dependence on A drops 
out and G turns into the Dirichlet propagator 

Gxyix, x') = Gb(a; - a;') -G^{x-x' ^ 2z'n) . (66) 

Both properties, the Dirichlet boundary condition as 

well as the A independence for c ^ 00, carry over to 
each individual Feynman integral of the bare (dimension- 
ally regularized) correlation functions G^^'^\ Thus they 
hold at least to any finite order of perturbation theory. 

Now, suppose we compute the renormalized functions 
Gren'*^^ usiug RG-improved perturbation theory. The 
RG maps these functions with given initial values u, c, 
and A to their RG images at points on the RG trajec- 
tory with the values u{£), c{£), and A(^). To obtain the 
asymptotic behavior of the G^m'^'' in the large length- 
scale limit £ 0, we may replace the running inter ac- 
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tions constants u, c, and A by their hmiting values u* , 
'^md — ^^'^ for ^ provided the RG images of 
the correlation functions remain finite and nonzero in this 
limit. Below the upper critical dimension d*{m) (i.e., for 
e > 0), there is no reason to expect any problems in the 



limits 



and A ^ Al: According to RG-improved 



perturbation theory, neither A — A!!j_ nor u — u* should be 
dangerous irrelevant variables. ^^'^^ 

On the other hand, as c 
with M > should vanish. To see this, note that the 
boundary condition (4) yields 



oo, all functions G^^n^^ 



0® = Cp ^ a„0p 



(67) 



for the Fourier transform of the bare boundary operator 
(/>®, where Cp was defined in Eq. (8). Using this and 
expanding in powers of 1/c, one sees that the leading 
contribution of the bare cumulants G^^'^^^ in the limit 
c —+ oo (with all other interaction constants kept fixed), 
is given by 

GWAf)(^.^^ = '^A/,2 5v,oc-i [l-(A/c)5„aJ(5(ri2) 



+ c-^^G(^'^)(a;;r), 



(68) 



where ri2 denotes the displacement ri — r2 of the two 
surface points in the case {N, M) = (0, 2). The functions 
on the right-hand side are defined through 

, iV M 
= l k—1 ^ 0=00, A=0 

(69) 

The extra term for {N, M) = (0, 2) in Eq. (68) is pro- 
duced by the zero-loop contribution to G^^'^^. Equation 
(68) generalizes the m = result (3.133) of Ref. 28 to 
the m > case. 

The large-c behavior (68) of the bare functions G^^''^^ 
implies that the RG images of their renormalized coun- 
terparts Gi^n^^ vary ~ cJ^(-i+o(e)) as c ^ oo, aside 
from extra terms proportional to S{ri2) and dada S{ri2) 
for {N,M) = (0,2). All in all, the following conclusions 
can be drawn: 

(i) The A dependence that remains when corrections 
to scaling ~ A — A!j_ are ignored (upon making the re- 
placement A ^ A^) is restricted to nonuniversal ampli- 
tudes (which are expressible via appropriate RG trajec- 
tory intcgrals^'^). 

(ii) The renormalized cumulants G^en^^'' approach zero 
in the limit c — > oo whenever M > 0, and satisfy Dirichlet 
boundary conditions. In other words, taken at c ^ oo, 
they vanish if one or several positions Xj of the operators 
<pren{xj) approach the surface.^^ 

(iii) The critical behavior of the G'^'^'^'^ at the LP or- 
dinary transition must match that of the G^''^\ and 
hence should be derivable in much the same way as in 
the m = case by a direct RG analysis of the latter 
functions. 



B. Renormalization of the cumulants G^'*'^ 

To set up such an analysis, let us consider the renor- 
malization of the cumulants G^'^\ Adding source 
terms to the Hamiltonian, we introduce the action 

A[4);J,Ji,oo] = dVJ(t)+ dAJi^^-dncj) 

- W[0]|,^^,^^o . (70) 

in which the Hamiltonian is given by Eqs. (l)-(3) with 
c = oo and A = 0. As before, we need the bulk 
reparametrizations (16). The surface operator 9„</> is 
multiplicatively renormalizable; we introduce its renor- 
malized counterpart (9„0)ren via 



1 /2 

a„0 = [^^(u, e) Zi,oo(u, e)] {dn<p) 



(71) 



Since the Fourier transform &'^''^\p) has /U-dimension 
one, it may contain additional primitive uv divergences 
localized on the surface of the form /i (Co + pCi) + 
c^^^ PaPaG2, where the renormalization functions Co, 
Ci and C2 would diverge as Co ~ A, Ci In A, and 
C2 ^ In A if we regularized the theory's uv singulari- 
ties by a large-momentum cutoff A. Thus the generating 
functional of the renormalized functions G^i-en^ can be 
written as 

^cx),reii[*^ , «^l,oo] 



= lnTr<^exp|>t 4>; J / ^/Z^, Ji^oo/ ^/Z^Z^c 



+ / dA 



M(Co + pCi)j2^ 
ic2ai/2^(a„Ji,^)2 I 

a=l -I J 



(72) 



where Tr^ means functional integration over 0. Hence 

we have 

G^rfn^ = [/x(Co + pCi)-ai/2C2a„a„]4'°5(ri2) 

^^_(jV+M)/2^-M/2^CAr,M) (73) 

for the renormalized cumulants in position space. 

In our dimensional regularization scheme, Co is finite 
and hence can be chosen to vanish. The In A singularities 
of Ci and C2 translate into poles at e = 0. For conve- 
nience, we fix Zi^oo and C2 by requiring that the poles of 
G(i'i)(p,z)|p=o and the additional ones of G^^'^'' {p)\p=o 
be minimally subtracted, respectively. Likewise, Ci 
could be determined by requiring that the additional 
poles of [dp&^'^\p)]p=o be minimally subtracted. How- 
ever, the explicit two-loop expression for Ci will not be 
needed in the sequel, so we refrain from determining it. 

The graphs of Gto'^'(p, z) to two- loop order are shown 
in Fig. 3. To this order, the only graph contributing to 
C2 is the analog of the last one, namely ®-|— rj-VH-® • In 
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G'oo'^\p,z) = o- h® + o~--^--|-® + 0-^-'='-'^-'^-f® 

/"\ 

FIG. 3: Graphs of G^'^\p, z) to two-loop order. The dashed 

line denotes the free Dirichlet propagator Gd- The vertical 
bars signify normal derivatives dn- Open and crossed circles 
indicate external points off and on the surface, respectively. 



Appendixes A and B the Laurent expansions of these 
graphs for /9 = r = are computed to the required order 
in e. Using these one arrives at the results 



1 



n + 2 u [11. + 2)[i) + 5) 



+ ■ 



3 2e 
n + 2 ,. 



36 



3 [ii(m)-J„(m)]-+0(«3) (74) 



and 



n + 2 



2 J2,3(m) 



3 6-m 



ml 



+ 0{v^). (75) 



Here bm is the constant introduced by Eq. (B3), 
/^^^ (m) and 7^^' (m) are the integrals defined by Eqs. (B5) 
and (B6), respectively, while J2,3(m) denotes a particular 
one of the integrals (49). 

The quantity 

ji (to) =Bm dy y""-^ ^m,d' iv) / dv 4>^_rf. (v) 
Jo Jo 

(76) 

is a new integral, not encountered in the bulk case. It is 
not difBcult to see that it can be reduced to a single inte- 
gral. Rewriting the double integral (76) as J^dv J^dy . . . 
gives 



jl (to) =Bra dv rf. (u) ilm,d' (^^) 



with 



9., 



J V 



dyy^'-^'^mAy) ■ 



(77) 



(78) 



The latter integration can be performed analytically; at 
d = d* (to) the result becomes 



8(TO-2)r(i + f) 

(m-4)r(2i±2) 

2(3Tn/2)-7 p( m-2 ) 



sin(TO7r/4) 



(79) 



If TO = 2 or TO = 6, then even the u-integration remain- 
ing in Eq. (77) can be performed to determine ji{m) an- 
alytically because the scaling function ^m.d* reduces to 
an elementary function. ^^'^^ Analogous statements ap- 
ply to the integrals J2,3(to), I^^\ni), and /(^^(to) (cf. 
Appendix E). One obtains 



ii(2) = i 



In 3 



ii(6) 



2 „, 27 
-- + 2\n — , 

3 16 ' 



. 2 J2,3(2) 

3 6 - TO 



m=2 



1 

54 



, 2J2,3(6) 



3 6-TO 



m^6 



91n(4/3)-2 
27 



(80) 



(81) 



(82) 



(83) 



b2l^^\2) = ^[Li2(2/3)-Li2(l/4)-ln(4/3)ln2] 



TT^ 1 

+ - ln(3/2) , 

192 12 ^ ' ' ' 



ln(4/3) 



m— *-6 



(84) 



(85) 



62/(2) (2) 
and 



TT^ + ln(4/3) [3 ln(4/3) - 8] 6Li2(l/4) 



b„,I^^\m) 



96 



1 - ln(4/3) 



m^6 



(86) 



(87) 



where it should be noted that bm/ (6 — to) approaches the 
finite value 2^^ 7r^/3 as to — > 6. 
The limit 



ji{0+)= lim ji(TO) = i 

m— ►0-1- Z 



(88) 



can be determined in a similar manner as the to ^ lim- 
its of the integrals j^, . . . , ju were in Ref. 15. For other 
values of to, no analytical results for ji (to) and the other 
integrals are available. We therefore resorted to numeri- 
cal means, proceeding in a manner analogous to the one 
described in appendix E of Ref. 15 (see Appendix E be- 
low). The resulting numerical values of ji{m), J2,3(to), 
I^^\m), and I^^^{m) are gathered in Table I, along with 
those for the (previously computed^^) integrals j^{ni), 
ja{m), jp{m), ju{m), and J„(to). 
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TABLE I: Numerical values of the integrals j,p{m), ja-im), jp{m), ju{m), Ju{m), ji{m), 100 x | J2,3{'rii), 10 x 6^/^^'(m), 
and 10 X bml'-^\m) 



m 


1 


2 


3 


4 


5 


6 


U (»") 


1.642(9) 


1.33333 


1.055(6) 


0.803(7) 


0.57(4) 


0.36521 


jcr (to) 


1.339(4) 


4.74074 


10.804(3) 


20.067(7) 


32.95(4) 


49.77778 


ip(m) 


0.190(6) 


0.88889 


1.999(9) 


3.464(1) 


5.23(4) 


7.26958 


ju (to) 


-0.203(7) 


-0.40547 


-0.624(2) 


-0.880(1) 


-1.21(1) 


-1.72286 


Ju(to) 


0.383(8) 


0.28768 


0.200(8) 


0.119(8) 


0.04(3) 


-0.02971 


ji (to) 


0.47289 


0.45069 


0.43092 


0.41273 


0.39577 


0.37983 


3 6-,r •^2,3(to) 


1.66642 


1.85185 


1.98142 


2.07349 


2.13819 


2.18198 


inf),„ /<i'(ni) 


n.:ifi44(-'')) 


0.40504 


n,4:i:!7(l) 


0.454:i8 


0.4f)90(n) 


n. 47047 




(1.8,S2(2) 


().!),S1.-)U 


1,0. Hi) 


1.110(7) 




1,18720 



C. RG equations and scaling 

Exploiting the relation (73) between the bare and 
renormalized cumulants in a standard fashion yields the 
RG equations 



Woo) ^ 



{N + M)r]^ M7?i,c 



rv^ von 



2 2 
= 4'° (M7^l - 7^2 d^d^)S{r,2) (89) 



with 



(90) 



7^l(/z,p,w) = Z^^Z^^^ 3^lo[-^<?i-^i,oo/u(Co + pCi)] , 

(91) 

and 

n2{u,u) = a-^I^Z^^Z^)^ iidX[Z4>Zi,oo<y^'^C2] , 

(92) 

where r?i,oo(^^) is defined by Eq. (29) with p = (l,oo). 

The RG equation for G^' ren is inhomogeneous because 
of the additive countcrtcrm in Eq. (72) and the imphed 
additive term in the reparametrization relation (73). If 
Co, Ci and C2 are determined by minimal subtraction of 
poles (so that Ci and C2 are Laurent scries in e with a 
vanishing regular part while Co = 0), then the uv finite- 
ness of TZi and 7^2 together with the familiar expressions 
(33) and (35) for and the RG functions 7]^ imply that 
Eqs. (91) and (92) simphfy to 

TZi = —p udu Res Ci (w, e) , 1^2 = —udu Res C2 {u,e) . 

e=0 e=0 

(93) 

Thus, to order , TZ2 is simply given by — 2e times 
the term oc u^/e on the right-hand side of Eq. (75). For 
the exponent function ?7i^oo(m) we find from Eq. (74) the 
result 

71+2 

Vi,oo{u) = ^ u {l + u[ji{m) - Ju{m)] } + 0{u^) . 

(94) 

We arc now ready to exploit the RG equations (89) 
by an appropriate generalization of the analysis given in 



Ref. 14. Owing to the p dependent term of the func- 
tion (3t- [ci. the set of equations (33)], the flow equations 
(40) for f and p are coupled. We therefore introduce a 
nonlinear scaling field 

9t{t,p,u) = T + c^2 (u) (95) 

such that the associated running variable is a solution to 

e^^gr{^) = -['^ + Vr{u)]gr{i), 9r{^) = 9r ■ (96) 

Substituting the ansatz (95) into the flow equations (40) 
leads to the condition 

f3uiu,e) + br (u) = [77, (u) -2inp(u)]cl.{u) . (97) 

which is easily solved for u near u* to obtain 



= ^^^^ 

Since we shall ignore corrections to scaling u — u* below, 
it will be sufficient to keep the flrst term on the right- 
hand side. Our one-loop result (39) for 6t(u) yields for 
it: 

c^(m*) = = --+Oiu*). (99) 

Let us define the RG trajectory integrals 



[u, u] = exp • 



(100) 



L = (j), a, p, T, (1,00), 

and the nonuniversal constants 

E*{u) = E,{u*,u) , i = 4>, a, p,T, {1,00), (101) 

they approach in the infrared limit £ ^ 0. Then the 
solutions to the flow equations (40) and (96) for a, p, 
and Qr can be written as 

a{i) = r'''- E4u{e),u]a ^^J-'''' E;{u)a , (102) 
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and 



p{£) = £-(-^+ip) Ep[u{£),u] p ^« J"^/"'^ E*{u) p , (103) In the special case of G£-,?L(ri2), the inhomogencity 

implies a contribution (x 6{ri2) linear in p we suppressed 

in Eq. (Ill), regarding G^'^lniri2) as a conventional 
function (rather than a distribution) with ri2 7^ 0. In 
momentum space this term yields a regular contribution. 
It reflects the fact that the momentum-dependent surface 
susceptibility xii{p) [cf- Eq- (52)] does not diverge at the 
ordinary transition and may be viewed as analog of the 
term oc ^1- (68)- 

Owing to the relations (68) between the bare cumu- 



(104) 

where we have introduced the familiar correlation-length 
exponent 



ui2 = {2 + n;)-' 

and the crossover exponent 

ip = vi2 (1 + r?*) . 



(105) 



(106) 



We choose £ = £r such that Qri^r) = ±1 for > 0, 
consider the limit Qr 0±, and introduce the scaling 
lengths^^ 



^l2 = ^X-'/ir^ [E;{u)\gr 



(107) 



and 



64 = 



-,1/4 



(108) 



where vn denotes the correlation-length exponent 

2 + 7?; 



= 61/12 = 



4(2 + 7?*) 



(109) 



The RG equations (89) can now be solved in a straight- 
forward manner via characteristics. It is convenient to 
utilize instead of r the nonlinear scaling field Qt, consid- 
ering the cumulants G^ren'' as functions of p, gr, and u. 
Upon substituting the trajectory integrals (100) by their 
fixed-point values (101) and the running interaction con- 
stants CT, p, and Qt by their asymptotic expressions dis- 
played on the far right of Eqs. (102) (104), one finds that 

the Gr^'*^' behave near the Lifshitz point as 



{ra,ri3,z; p, gr, u, a, p) 



, {N+M)/2 



-2-')i*,o 



EL 



■^12 



M/2 



?;4 C,12 C,12 



(110) 



Here and stand for the sets of all a and /3 coordi- 
nates parallel to the surface, respectively, while z repre- 
sents the set of all coordinates perpendicular to it. Fur- 
ther, G^±^^ is a universal scaling function, given by 

G^r^,'f\rc,rfi,z;p) 
= G^^;J^i^[ro,,ri3,z;p,gr=l,u*,l,l]. (Ill) 



lants G^^'^^ in the limit c 00 and G^''^'' the scaling 
forms (111) carry over to the behavior of the former at 
the LP ordinary transition. An immediate consequence 
is that the critical exponents at the LP ordinary transi- 
tion can be expressed in terms of four independent bulk 
critical indices, e.g., the correlation exponent 



Vl2 = % 



(112) 



and the previously introduced anisotropy, correlation- 
length, and crossover exponents 6, V12, and (p [cf. 
Eqs. (105) and (106)], and a single surface critical in- 
dex. As the latter one can choose the critical exponent 
^ord governing the asymptotic dependence 



mi(T,p = 0) 



(113) 



of the surface magnetization mi at p = (and given 

values of the remaining interaction constants c, u, and 

A). From the behavior of G^.i-en predicted by Eq. (Ill) 
one reads off that 



P^'^ = m [d + m 



m{e-l)]/2. (114) 



D. BoundEiry operator expansion 

In the above derivation of the scaling behavior of the 
cumulants G^^'^^ at the ordinary transition we relied on 
the large-c form (68) of these fimctions. An alternative 
way to arrive at the same conclusions, which gives further 
insight, is the use of the boundary operator expansion 
(BOE).^*''^^ Since both the bare and renormalized theo- 
ries satisfy Dirichlet boundary conditions for c = 00 and 
c = 00, there is no term cx (f)^ contributing to the BOE 
of <p{r, z). Just as in the CP case, the leading operator 
should be 9„0. Hence we write 



(115) 



Consistency with the RG equations (89) requires that the 
function Goo (2) must satisfy the RG equation 



(p(.°°)-r?i,oo/2)Gco(z) = 0. 



(116) 



Solving this (with p set to zero) yields a short-distance 
singularity for z <^ ^12 of the form 



Goo(2) ~ const (a<2)1+<~/^ 



(117) 
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By analogy to the CP case, the exponent is governed by 
the difference between 



I3i = vi2 [d-2 + Tji2+m{e-l)]/2, 



(118) 



the critical index of the bulk order parameter mb = 
{(t>\z=oo), and /35"'<^; we have 

{PT^-l3i)/vi2 = l + r^U2. (119) 



E. Exact critical exponent of the surface energy 
density 

In the foregoing subsections IV A-IV D we have shown 
that the critical exponents of the cumulants G'^^'^^ at 
the ordinary transitions can all be expressed in terms of 
four bulk critical indices (e.g., 77(2, 9, 1^12, </?) and a single 
surface critical index such as /JJ'^. We now turn to a 
discussion of the behavior of the surface energy density 



£:i^<[0^(r)]V2) 



(120) 



at the LP ordinary transition. 

The surface energy density at the CP ordinary transi- 
tion is a well-known example of a local boundary density 

whose critical exponents are expressible in terms of bulk 
exponents: Its leading thermal singularity is of the same 
form |r|^~" as the bulk free energy. ^^^^ Our aim here is 
to show that the same holds true for the surface energy 
density at the LP ordinary transition. 

A particularly c;asy way to obtain this result is to argue 
along the lines of Rcf. 49 and consider a finite system that 
extends in the z-dircction from z = —L2 io z — Li, has 
cross-sectional area A, and different values fi, pi and T2, 
P2 in the regions 2Ji with z > Q and 2J2 with 2; < 0, 
respectively. Its Hamiltonian reads 



H2 



xCb{x;fi,pi) 



d'^xCb{x;T2,p2) 



2J2 



(121) 

where Ch{x; f, p) is the bulk density (2). We now choose 

values of T2 and p2 such that the bulk at large — z > is 
not critical and in the disordered phase when the thermo- 
dynamic limit is taken. ^° Since the half-space z < for 
the chosen values of f2 and p2 is neither critical nor has 
no long-range order, its net effect is to generate modified 
effective short-ranged 0{n) invariant interactions at the 
interface z = 0. The boundary critical behavior which 
occurs at this interface as fi — > tlp and pi — > plp should 
be in the universality class of the LP ordinary transition. 

Suppose the interface is displaced a small distance AL 
downwards, so that the heights Li and L2 of the re- 
gions 2Ji and 9J2 increase and decrease by AL, respec- 
tively. The change of total free energy per volume A AL 
caused by this displacement in the limit A 00 and 
Li, L2 ^ 00 becomes 



/b(Tl,Pl) - /b('r2,P2) 



= (n - h) Si + 



Pi - P2 



Y: ({d.ct>'^r) . (122) 



Here the left-hand side is simply the change of bulk free 
energy per unit volume caused by the alteration of values 
{t2,P2) {ti,Pi) within the region — AL < 2; < 0. The 
right-hand side follows by expressing this change in terms 
of 



AH 



i-AL t- 
J ^^J^^ 



Tl - T2 



0^ + 



Pi - P2 



2 " 2 ^ 

a — 1 

(123) 

the corresponding change of the Hamiltonian, as 
-(AAL)-iln(e-^'«) « {AH)/{AAL), where (.) indi- 
cates a thermal average with the Boltzmann factor e^^^ 
of the original unperturbed system. For small AL we 
have J^^dz f{z) w /(O) AL. Assuming translational in- 
variance parallel to the interface then yields the right- 
hand side of Eq. (122). (This assumption excludes sit- 
uations in which the order parameter varies along the 
interface. In such cases of broken translational invari- 
ance along the interface the right-hand side of the identity 
(122) must be replaced by its spatial average A~^ J dA. . . 
over the interface.) 

For small deviations fi — tlp and pi — pLp from the 
LP, i.e.. small values of the scaling fields gr [cf. Eq. (95)] 
and p, the first term on the right-hand side of Eq. (122) 
behaves as 



/b(Tl,Pl) « \gr 



reg . 



(124) 



Here "reg" stands for regular contributions, and the plus 
and minus signs of the subscript ± refer to the cases 
5r > and 5t < respectively. 

Noting that ^(0) 0, we can equate the leading sin- 
gular parts of the identity (122) to obtain the asymptotic 
behavior 

o o m 

I o Nosing , PLP — P2 sr^ 1,0, .<8N2\sing 

(tlp - T2) f 1 > ((5„0 ) ) 

~^±(0)|t|^-"' (125) 

for T — > 0± at p = 0. Hence one of the two surface quan- 
tities on the left-hand side must contain the singularity 
on the right-hand side. 

The critical line fc = /cl (p) is not a straight line in the 
vicinity of the LP. On the section separating the disor- 
dered phase from the modulated ordered one, its asymp- 
totic behavior near the LP is governed by the crossover 
exponent i.e., fc — tlp ^ [p — pi^'pY^'^ ■ Nevertheless, 
choices of the point {J21P2) in the disordered phase with 
P2 = Plp and T2 ^ tlp should be possible. Since for such 
choices only the term proportional to the surface energy 
density Si contributes on the left-hand side of Eq. (125), 
it is clear that Si must have the bulk free-energy singu- 
larity on the right-hand side. 

This result means that the scaling dimension A[(0®)^] 
of the surface energy density operator (0®)^ coincides 
with that of the volume / dV, i.e.. 



A[(</.'^)^] =d + 0(m-l) . 



(126) 
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As usual, the amplitudes ^±(0) are nonuniversal, but 
the ratio J^+(0)/J^_(0) is a universal bulk amplitude ra- 
tio. The corresponding amplitude ratio of £i must have 
the same value. 

There is an instructive alternative way to confirm this 
result. Knowing that the ordinary fixed point T'*^.^ is 
located at c = oo, one can study the effect of the Hamil- 
tonian's boundary term (c/2) tfp for large but finite c 
to leading order in 1/c. 

Expanding the free propagator (7) in powers of 1/cp 
gives 

G{p\z,z') - Goip;z,z') 
= c-^e-^^^'-''^ +0{c-^) (127) 
= (Gd5„)(p; z, 0) (a„GD)(p; 0, z') + 0{c-^) , 

where i9„ acts to the left. More generally, one notices 
that the expansion of the free propagator in powers of 
l/cp in the present m ^ case is completely analogous 
to the expansion of its m = counterpart in powers of 
1/c, with identical expansion coefficients [cf. Eq. (11) of 
Ref. 47]. 

If we restrict ourselves to first order in c~^, we can 
make the replacement c"^ — > c~^. In analogy with the 
m = case the expansion of the iV-point cumulant to 
first order in c^^ therefore becomes 

where G^'^'^"^ means a cumulant involving N fields 
(j>a and a single insertion of the boundary operator 
J^dA (9„0)^/2, evaluated at c = oo and A = 0. Proceed- 
ing as in Appendix C of Ref. 47, one can easily derive the 
relation 

AT 
i=l 

(129) 

where Zi is the ^;-component of Xi. Since the renormaliza- 

tion factor Z^^^^ of the cumulant on the left-hand side 
must also renormalize the one on the right-hand side, 
the inserted boundary operator needs no rcnormaliza- 
tion. Thus the associated interaction constant ^ 
must scale naively as l/ji at the ordinary fixed point. 
Noting that the RG oigcnexponcnt jji of a boundary scal- 
ing field gi and the scaling dimension A[C'i] of a local 
boundary operator 0\ (in the case of the parallel surface 
orientation considered in this paper) are related via 

A[0^]=d-l + m{e-l)-yi , (130) 

one immediately concludes that the scaling dimension of 

A[(5„0)^] of the inserted operator (9„0)^/2, which rep- 
resents the surface energy density at the ordinary transi- 
tion, is indeed given by (126). 

The above findings indicate that the boundary opera- 
tor with smallest scaling dimension contributing to the 



BOE of the local energy density operator (t>'^{x) (be- 
sides the one-operator 1) is T®, the zz-component of 
the stress-energy tensor T^^, taken at the boundary (cf., 
e.g., Ref. 51). 

V. CRITICAL EXPONENTS OF THE 
ORDINARY TRANSITION 

A. Definition of surface correlation exponents 

In order to characterize the power-law decay of the 
pair correlation function G*^^'"^ (xi, 032) of semi-infinite 
systems at a critical point in the limits where the separa- 
tion X1 — X2 increases parallel or perpendicular to the sur- 
face, one conventionally introduces analogs of the usual 
correlation exponent r} termed r/y and tjj,, respectively. 
In the case of anisotropic scale invariant systems with 
boundaries we must be more careful because distances 
along different axes scale differently and the orientation 
of the surface matters. 

Clearly, distinct sets of surface correlation exponents 
'ni2,\\i '^l2,^_, and f]i4,||i ^i4.±, could be introduced in anal- 
ogy with the bulk correlation exponents 77(2 and rju. How- 
ever, this is unnecessary since the "Z4 exponents" are sim- 
ply related to their "Z2 counterparts" via the anisotropy 
exponent 9, just as the introduction of rju could be 
avoided by writing the exponent 4 — r^n as 6 {2 — r]i2). In 
order to define the exponents ?7|| = 77(2, || and r]± = 77/2, _Lj 
consider the pair correlation function G^^'°)(r; 2, 2') at 
the Lifshitz point between two points x = {r,z) and 
x' = (0, z'), and let us introduce 

f = y/raTa and R = sjr^ , (131) 

the lengths of the components (r^) and (r/j) of their par- 
allel separation r. In the limits ^; — > 00, i? — > 00 or 
f — > 00, with z' and the respective other distances kept 
fixed, this function decays as 

R-ld-2+m{e-l)+r„] [ R ^ ^o' , 
._[d_2+m(e-l)-Fr)||]/0 ^ f.^QQ _ 

(132) 

These relations translate into the small-momentum be- 
havior 

^ . ^si„g j P'"-' for (Pa) = (0), ,.„„x 

-jpC'/ii-i)/^ for (p^) = (0), ^^^^> 

and52 

V /P""-' for(pj = (0), 

of the local surface susceptibility (52) and the layer sur- 
face susceptibility 

/■oo 

Xi(p)= / G^'''\p;z)dz, (135) 
Jo 
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respectively. 

The large-c behavior (68) in conjunction with the scal- 
ing forms (111) implies that 



2 + Viz + V*i, 



(136) 



The corresponding result for tj'^^ is equivalent to the scal- 
ing law 



V± = {vi2 + V\\)/'2 



(137) 



which is well known from the CP case.^®'^^'^^'^^ 

Let us also note that the decay law of G^'^'^^ (involving 
1]'^'^) displayed in the first line on the right-hand side of 
Eq. (132) carries over to the case where x moves away 
from the surface along an arbitrary direction not parallel 
to the surface. This can be shown as in the CP case,^*'^'^ 
utilizing the BOE (115) and Eq. (117). 



B. € expansion of surface critical exponents and 
estimates for three dimensions 

We now turn to the e expansions of the critical expo- 
nents at the ordinary transition. Substituting the e ex- 
pansion (42) of the nontrivial root u* of /?„ into Eq. (94) 
yields 



n + 2 ^ {n + 2) 



+ 2 



n + 8 + 
20-h 7n 



16 {m + 2) 2 

2 



+ 0(e3). 



Ju{m) + Q — —rji(m) 



n + 2- 



(138) 



For m — > 0, m = 2, and m = 6, the exactly known 
values^^ 



j^(0) = 2, i<,(0) = 0, J„(0) = l/2 



4 128 4 

i^(2)=3, >(2) = — , Jn(2) = ln- 



(139) 



and 



jm = 3 [1-3 In -J 
J«(6) = ^ - 31n^ , 



• 448 



(141) 



can be inserted into Eq. (138), along with those of ji 
given in Eqs. (88) and (80), to convert it into an analyti- 
cal expression. Doing so one finds that the result (138) re- 
duces in the limit m to its established analog"'*'"^^''''^ 
for the CP ordinary transition, Eq. (IV. 35) of Ref. 37, as 
it should. 

We can now replace 77/2 in Eq. (138) by its known O(e^) 
expression [given in Eq. (85) of Ref. 14 and Eq. (61) of 



Rcf. 15] to obtain the e expansion of 77°'''^ to order e^. 

The analogous "direct" series expansions to order of 
the surface exponents (3°^'^, vT'^^ and the surface suscep- 
tibility exponents 7^"''^ and ^'fl'^ follow in a similar fashion 
by combining the e-expansion results for the bulk expo- 
nents given in Eqs. (61)-(66) of Ref. 15 with ours for 77^ ^ 
and the scaling relations (137), 



(3i = ui2[d-2 + 9{m-l)+^]/2 



711 = 1^12 (1 - m) 



and 



71 = 1^12 (2 - V±) , 



(142) 



(143) 



(144) 



all of which arc straightforward consequences of the 
scaling forms (111), the large-c form of the cumulants 
G^^'^\ and the short-distance behavior (117). They 
can be derived in much the same way as in the m ~ 
for the parallel orientation of the surface con- 



case 



28,37 



sidered in this paper. 

In Table II wc have gathered the numerical estimates 
of various critical exponents of the LP ordinary transition 
one obtains for the uniaxial one-component case in three 
dimensions by cvahiating their e expansions to zcroth, 
first, and second order at e = 3/2. The critical indices 
/3s and 7s, defined through the behavior rUs \t\^' and 
Xs ~ kl"^" for r ^ 0, p = 0, of the excess magnetization 

POO 

ms= / dz[G^^'°\z)-G^^'°\oo)] (145) 
Jo 

and the excess susceptibility 

pOC POO 

Xs = J dz j dz' [G(2.o) (p.z,z')- (p; 00, 00)] 

(146) 

do not involve an independent surface exponent because 
they can be expressed in terms of bulk exponents as 



(140) and 



/3s = A - n2 



Is = 71 + n2 



(147) 



(148) 



We have included them in Table II to compare also them 
with the Monte Carlo results of Ref. 21. 

The direct series expansions to O(e^) of the exponents 
listed in Table II are reasonably well behaved. We there- 
fore expect the corresponding estimates to be fairly reli- 
able. Yet for some of the series, Pade [1,1] approximants 
give significantly different values. In other, more favor- 
able cases the differences are small. One example of this 
kind is j3°'"^ . The Pade [1,1] approximant and the solu- 
tion to the scaling relation (119) yield 



/?r^(3,i,i) 



0.683 , Pade [1, 1] , 

0.628 , with bulk exponents 



(149) 
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TABLE II: Estimated values of critical exponents at the LP 
ordinary transition in d = 3 dimensions for m = 1 and n = 1 . 
The values in the columns labeled "MFT", "0(e)", "0(e^)", 
and "MC" correspond to mean field theory, to the estimates 
obtained by setting e = 3/2 in the e expansions to first and 
second order, and to Monte Carlo results, respectively. 



MFT 0(e) O(e^) MC 





2 


1.5 


1.133 






1 


0.75 


0.586 






1 


0.75 


0.697 


0.687(5) 




1/2 


0.75 


0.947 


0.82(4) 




-1/2 


-0.375 


-0.212 


-0.29(6) 


/3b 





-0.375 


-0.462 


-0.46(3) 


7s 


3/2 


1.875 


2.106 


1.98(8) 



"Taken from Ref. 21. 



respectively. To obtain the second number, the "best 
estimates" f3i ~ 0.220 and 1^12 ~ 0.746 of Ref. 15 were 
inserted into Eq. (119) together with the direct series 
estimate V* 00 — —0.9060. Both so-obtained numbers 
are approximately within 10% of the estimate given in 
Table II. 

The exponent 7°'''^ provides an example of a case 
in which the = 3 value predicted by the Pade ap- 
proximant, namely 7°'''^ ~ 1.670, deviates considerably 
from the direct series estimate. However, the relation 
7°'''^ = (^(2(1 — V12 ~ Vi oc/2), applied with the just men- 
tioned values of 1^12 and rjl ^ and the "best estimate" 
r]i2 ~ 0.124 of Ref. 15, yields 7f ^ ~ 0.99, which is rea- 
sonably close to the direct series estimate (Table 11) . We 
therefore discard the Pade [1,1] estimate in this case. 

As can be seen from Table 11, our estimates are in 
good agreement with Pleimling's^^ Monte Carlo esti- 
mates. Normally the Monte Carlo results for Pi are more 
accurate than for susceptibilities such as 71 and 711. It 
is therefore gratifying that the agreement is particularly 
good in the case of 

Of interest is also the case of the ordinary transition 
at the uniaxial Lifshitz point of the axial next-nearest 
neighbor XY model (ANNNXY model), with n = 2 and 
m ~ 1. This model was recently argued to describe the 
bulk critical behavior of Tb.^^ To our knowledge, nei- 
ther experimental nor theoretical results for the surface 
critical exponents at the corresponding ordinary transi- 
tion have been published. On the other hand, Monte 
Carlo results for the bulk critical exponents of the three- 
dimensional ANNNXY model have been obtained a long 
time ago.^^'^'' As discussed in Ref. 15 (see its tables 4 
and 5), Selke's values'" A = 0.2±0.02 and 7; = 1.5±0.1 
are in reasonable agreement with the O(e^) estimates 
Pi ~ 0.276, and 7; ~ 1.495 obtained there. Inserting 
the latter together with the O(e^) result V12 ~ 0.757 into 
Eqs. (147) and (148) yields ~ -0.48 and 7^ ~ 2.25. 

To obtain estimates of p'^'^ for d = 3, m = 1, and 
n = 2, we exploited our above results for its series expan- 
sion to order . Direct extrapolation of this (truncated) 



scries and Pade [1,1] and [0,2] approximants yielded the 
values Pf^ ~ 0.72, 0.71, and 0.75, respectively. Aside 
from experimental tests of these predictions, checks via 
high-precision Monte Carlo simulations would be very 
welcome. 



VI. SUMMARY AND CONCLUSIONS 

We have investigated the critical behavior of semi- 
infinite d-dimensional systems at an m-axial bulk Lifshitz 
point. To this end we have considered systems whose 
bulk universality class is described by an n-component (j)^ 
Hamiltonian with the bulk density (2). Assuming that 
the surface plane of the system is oriented parallel to 
all potential modulation axes in the modulated ordered 
phase, we have constructed appropriate "minimal" semi- 
infinite extensions of these continuum models with the 
following properties: 

(a) They are compatible with the bulk universality 
class of Lifshitz critical behavior described by an infinite- 
space model whose Hamiltonian has the density (2). 

(b) They are minimal in the sense that redundant and 
irrelevant surface terms compatible with the presumed 
0{n) symmetry and short-range nature of the interac- 
tions and their surface-induced perturbations have been 
dropped. 

(c) Taken at their respective RG fixed points, these 

models represent the (surface) universality classes of the 
associated ordinary, special, and extraordinary transi- 
tions. 

The resulting contributions to the Hamiltonian (1) lo- 
calized on the surface correspond to the surface density 
(3). They differ from the usual one known from the 
much studied (m = 0) case of a critical point, namely 
(c/2)^^, by a similar derivative term. We have shown 
that the latter is required for rcnormalizability of the 
model. Physically, this means that such a term, if ini- 
tially absent, gets generated as short wave-length degrees 
of freedom are integrated out. The fluctuating boundary 
conditions implied by the Hamiltonian (l)-(3) are given 
in Eq. (4); in the pz representation they correspond to 
Robin boundary conditions with a momentum-dependent 
Cp [see Eqs. (4), (8) and (67)]. 

The boundary term oc A gives rise to a dimensionless 
interaction constant A, on which the surface coimterterms 
and renormalization functions depend for general values 
of c. We have clarified the fixed-point structure of the 
model, identifying the fixed points 7-**^^, 7^*^, and V^^ 
describing the ordinary, special, and extraordinary tran- 
sitions (cf. Fig 2), which are located at a nontrivial value 
(62) of A. 

In order to investigate the LP special transition, deal- 
ing with the A dependence cannot be avoided. However, 
as we have shown, this is possible in the case of the or- 
dinary transition by taking the limit c — > 00 from the 
outset. Just as in the CP (to = 0) case, both the regu- 
larized bare and the renormalized theories obey Dirichlet 
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boundary conditions if c = oo, and the single required 
surface critical exponent follows from the scaling di- 
mension of dn4>- Likewise, the proofs'* ^"''^ for the CP case 
showing that the surface energy density at the ordinary 
transition has a leading singularity of the bulk free-energy 
form ~ |tP~" carry over to the present m ^ one. 

Performing a two-loop RG analysis of the LP ordinary 
transition for general values of m, we have been able 
to determine the e expansions of the surface critical in- 
dices of this transition to second order in e. Extrapola- 
tions of these series expansions to = 3 dimensions yield 
values of the surface critical exponents for the uniaxial 
one-component case to = n = 1 in good agreement with 
recent Monte Carlo results^^ for the ANNNI model. 

There are several obvious directions for extensions of 
the present work. First of all, the special transition 
should be analyzed, preferably by means of a two-loop 
RG analysis giving the e expansions of its surface critical 
exponents to second order in e. Prom the results pre- 
sented above, the e expansions of the two independent 
surface critical exponents of this transition, the surface 
correlation exponent r/jj'' and the surface crossover expo- 
nent can be derived to first order in e with moderate 
effort. However, a detailed investigation of the special 
transition requires more work, is beyond the scope of the 
present paper and left to a planned subsequent one. Suf- 
fice it here to say that for our 0{n) symmetric Hamilto- 
nian with n > 1, the existence of a surface-ordered phase 
at temperatures T > is ruled out for bulk dimensions 
d < c??''"'' -|-1 = + m/2. Extrapolations of the e expan- 
sions of the special transition's critical exponents to d = 3 
thus make sense only in the scalar case, n = 1. For the 
interesting case m = n = 1, Monte Carlo results for the 
exponent (3l^ of the ANNNI model are available^* which 
may be used to check the resulting d = 3 extrapolation 
values. 

A second important line of extension is to consider a 
surface orientation perpendicular to one of the potential 



modulation axes (the "perpendicular case" of the Intro- 
duction). This case is somewhat more difficult to handle 
than the parallel surface orientation treated in this paper. 
Since the z-direction then scales naively as rather 
than as 1 //x, more potentially dangerous boundary terms 
(whose coupling constants have nonnegative RG eigen- 
exponents at the Gaussian fixed point) exist. Since the 
Landau-theory equations for the order parameter pro- 
file and the free propagator involve the fourth derivative 
d'^/dz^, two boundary conditions rather than the sin- 
gle one, Eq. (4) or Eq. (6), are required and found. ^"^ A 
proper study of this case beyond Landau theory*^ also 
requires a careful identification of redundant surface op- 
erators. 

Another obvious limitation of the present work is its 
lack of detailed investigations of the homogeneous and 
modulated ordered phases, the extraordinary transitions, 
as well as of profiles of the order parameter and other 
densities at the LP. The progress made here indicates 
that such studies, though technically demanding, should 
be possible, involving no fundamental new difficulties. 

Unfortunately, there are, to our knowledge, yet no ex- 
perimental results on surface c;ritic;al behavior at bulk 
Lifshitz points available with which we could compare. 
Owing to the richness of the physics involved, detailed 
experimental studies would be very interesting. We hope 
that the progress made recently on the theoretical side 
will foster such experimental activities. 
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APPENDIX A: LAURENT EXPANSION OF THE FEYNMAN INTEGRALS OF G! 



In this appendix, we briefly describe the calculation of the Laurent expansion of the one and two-loop graphs of 

the function G^x>'^\x,x') to the required orders in e~^. 

The Dirichlet propagator (66) to which the free propagator (7) reduces for c — oo can be written as 



Gb{p;zi,Z2) = 



-kp\zi-Z2\ 



-kp{zi+Z2) 



Gh{P; Zl - Z2) - Gs{p; Zi,Z2) 



(Al) 



in the pz representation. The part depending on |zi — Z2I is the bulk propagator Gh- The subtracted piece 
Gs{p; Zl, Z2) = Gb(p; zi + Z2) is the image term which ensures that the Dirichlet condition holds. 

As mentioned in the caption of Fig. 3, we use open and crossed circles to indicate external points off and on the 
surface 5B, respectively, a broken line to represent Gd, and perpendicular strokes to mark ^-derivatives dz. Accordingly, 
Eq. (Al) can be depicted as 



-O . 



O- -O = O— ^ — O 

Zl Z2 Zl Z2 Zl Z2 



(A2) 



The Feynman graph expansion of the two-point function G)xi {p\zi) to two-loop order is shown in Fig. 3. The 
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zero-loop term 



O-- 

Zl 



-h® = (Gd5„)(p;^i) = d,,Gj,{p;zi,Z2 = 0) = e'^-^^ 



(A3) 



follows directly from Eq. (Al). 

In order to determine the one and two-loop terms, we amputate the external free legs of the graphs and compute 
the Laurent expansion of the resulting amputated graphs, considered as generalized functions (distributions) of the 
dimensionless variable /xz and iJ.z' in the pz representation.^^ The amputated legs play the role of test functions^® on 
which the latter distributions act. Once the Laurent expansions of the amputated graphs are known, those of the full 
graphs can be computed in a straightforward fashion by performing the remaining integrations over z and z'. For 
convenience, we will set /z = cr = 1 in most of the sequel. Whenever necessary, the dependences on n and a can easily 
be reintroduced by dimensional arguments. 

The one-loop term of G^'^'(p; zi) reads 



0----'--HhS = -- 



n + 2 u 



dzGT,{p]zi,z) {Gudn){p;z)GD{x,x) , 



zi 3 2 

where the internal loop, Gr,{x,x), according to Eq. (13) is given by 

rt^nr.\- r(Oy^\- P r(2 - e) sin(e 7r/2) ^ 



(A4) 



(A5) 



Considering z"^ ^ as a generalized function, we can calculate its action on a test function (p{z) and expand the result 
about e = to obtain^^-^^ 



f 

Jo 



dzz'^ <f{z) 



1 d^p{z) 
e dz 



z=0 



+ {zf,^{z))+0{€) 



where z_^'^ is a standard generalized function, defined as 



/ -^W{z)-m-e{i-z)z^'{Q)\. 

Jo z 



Equations (A4)-(A5) yield 



O---— --H-® = U — - - [1 + Ru{p, zi) e + 0{e^)] o- 
Zi o 4e zi 



-y<s 



with 



and 



Rr){p,zi) = Ce - I + (^z^^,GTiip;zi,z) {GD^n)ip■,z)^/o- 



-[<S> 



= Ei(-2KpZi)e2^''^i -ln(2Kp) , 



where Ei(a;) is the exponential-integral function. 
We next consider the two-loop diagram 



(A6) 
(A7) 

(AS) 

(A9) 

(AlO) 

(All) 



Zl 



'u(n + 2y 




6 


lo 



dz Gd (p; 21,2;) (Gd dn){p]z) 



■.j^dz'Gj,{x',x') l±-^Gu{p';z,z')Gu{-p';z,z') 



Performing the p' integration associated with the lower loop yields 

z' 



> \ 

2 X : 1 



r(l-e) sin(e7r/2) 



[\z-z'r'' + {z + z'y-'' - 2 [max(z, z')]'"' } > 



(A12) 



(A13) 
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where the muhiphcation by 2 on the left-hand side compensates the graph's Une-symmetry factor of 1/2, so that the 
result is precisely the integral over p' (last integral) in Eq. (A12). 
The subsequent z' integration in Eq. (A12) is straightforward, giving 



/•OO '•n 

2 / dz' Gt,{x' ,x')\ ) = - 
Jo 



r(2-e) sin(e7r/2) 



'jm/i 



1 / 1 ^r(e)r(e-l) 



1-e tl-e r(2e-l) 



(A14) 



Using Eqs. (A6) and (A7), one can easily compute the action of this distribution on the external lines. The Laurent 
expansion of the result becomes 



The diagram 



1 1 1 I 







) 


■-1 1 + ^ 




. 3 , 




4e2 






2) 


2 /'CO p 

/ 




6 




Jo Jo 



8e 



+ O(e0 



O- h® 

zi 



(A15) 



xGd(p; z, z') Gu{x, x) Gu{x', x') 
is quite similar. Its resulting Laurent expansion reads 



/n + 2 



4 V 3 



1 _^l + 2^^^ ■ 

56^ 8e 



O- h® . 

Zl 



(A16) 



(A17) 



In deriving this equation we used the fact that the following integral of the generalized function z '^^'^z' e '^'^ ^1 
with a test function if{z,z'^ has the Laurent expansion 



/•OO /"OO 

/ dz \ dz'ip{z,z') 
Jo Jo 

^ {d,d,,F){0,0) 1 
2e2 + e 



«-^(5,i^)(o,z')) + 



(afF)(0,0)-(a2F)(0,0) 



+ 0(6°) 



(A18) 



where F is defined by 

F{z, z') = [<f{z, z') + ip{z', z)] . (A19) 

In addition, Eqs. (A5) and (A6) were employed. 

The computation of the diagram O r_ -l — f® is more cumbersome. In the position-space representation it 

can be written as 

71^ n + 2 



Zl 



6 3 



where A is an integral operator defined through 

{il),A^) = [ d'^x [ d'^x'i}{x)Gl{x,x')ip{x') 



with 



Since 



tlj{x) = Gu{xi,x), <f{x') = {GDdn){x' ,r2) ■ 

i^P, A^) = (Va A^)^{{^ + ^p,Ay + ^lj])- ((^, A^)- (V', ^ V')}/2 , 



(A20) 
(A21) 

(A22) 

(A23) 



it is sufficient to calculate the diagonal element {(p,A(f) for a general test function ip. Upon substituting Eq. (Al) 
into Eq. (A21), one obtains 



(A24) 



20 



where Ij is defined through 

Ij = {ip,Aj>f )= [ d'^x [ d'^x' ip{x)[Gb{x-x')f~^ [G,{x,x')y if{x') . (A25) 

Let us calculate the terms on the right-hand side of Eq. (A24) one by one, starting with the "bulk" graph 

Io = J dzj dz'Jd'^-\Jd'^-\'^{r,z)(p{r',z')Gl{r-r',\z-z'\). (A26) 

Exploiting the scaling property of the free bulk propagator, 

Gb(r;^) = Gb[(r„),M;0] = 1^1^-2 Gb[(r„|^|-i/2),(r;3/|^|;l] , (A27) 

and making the changes of variables pa = \z — z'\~^^'^ (r^ — ra) and pfj = \z — z' [r'p — r p) , we can 

recast Eq. (A26) as 

Jo = jd''-\ j^dz £dz' \z - ^{r, z) Jd^-'p (p[(r„ + VW^\Pa), {r0 + \z - z'\ p0),z'] Gl{p; 1) . (A28) 

This integral has a pole at e = caused by the ultraviolet divergence on the line \z' — z\ = which is due to the factor 
\z- z' 



'P*^ ^. For our purposes here it is sufficient to compute just the singular part of the integral (A28). To this end, 
we expand its integrand in powers of ^y\z — z'\ pa and \z — z'\pi3. Subsequent termwise integration over z' then yields 

lo = -Jd'' rj dzip{r,z)^— ^ — + d0d0ip{r,z) 

i 



4em 

where we have introduced 



/2(1) fd''-^rip{r,0)dadMr,0) + O{e'>) , (A29) 



/fe(t)^|yd'^*-V(p«'Pa')'^'<^b(P;l)Gb[(Pai'/'),(p/3t);l]} ^, fc = 0,2,4, (A30) 

and utilized the identity 

dz dz' \z - zT-^'xiz) X{z') = ^^^_^y_ I dzx{z)^^x{z) + O(e0) , (A31) 

valid for A'^ e N and test functions x{z)- 

Next, we turn to the terms Ij with j > 0. The associated Feynman diagrams contain the image part Gs{x, x') = 
Gh{r — r';\z + z'\), which may induce additional ultraviolet singularities^^ at ^; = and z' = 0. To compute 

h= f d'^x ( d'^x' ip{x)Gl{x-x')G^{x,x')ip{x') , (A32) 

we transform from the integration variables {z, z' , r'} to new ones, {v, t, p}, defined by 

z = v^^^, z' = v^-^, r'^ = ra + Vvtpa, r'0=r0 + vtp/3 , (A33) 
and then use again the scaling property (A27) of the propagator to rewrite Eq. (A32) as 

//•oo fl 
d'^-^r dvv^'-^ dtf-'^g{r,v,t) (A34) 
Jo Jo 

with 

g{r,v,t) = ip[r,v{l + t)/2] Jd'^''^ p ip[{ra + pc,) , {r/s + vt p0),v {1 - t)/2] x 

xGl{p;l)Gk[{t'/'pa),{tp0);l] . (A35) 
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Upon performing the integration over t to obtain 



dtt y[r,v,t) = — ' ' 



/ ^[g{r,v,t)-g{r,v,0)]+O{e) , 
Jo r 



(A36) 



we can do the integration over v in Eq. (A34), using Eqs. (A6) and(A7). Straightforward calculations then yield the 
result 



/o(0) 
2e2 



^(r,0) <p'(r,0) + 



2e 



[/oW-/o(0)] 



+ 



^ ip{r, 0) dc^dMr, 0) l^dt hit) + ^ {vf,^\r, v/2)) I + 0{e°) . (A37) 



The Laurent expansions of the integrals h and I3 can be worked out in a similar fashion. They read 

72 = yd'^-v|l(p(r,0)^'(r,0)^"dt/o(i) 



(A38) 



and 



73 = yd''-v|^<p(r,0)^'(r,0) + ^¥'(r,0)a„5„^(r,0)|+O(e°). (A39) 

Equations (A29), (A37), (A38) and (A39) give us the diagonal elements {(fjAjf) for a general test function (p. 
Utilizing these results in conjunction with Eq. (A23), one can readily determine the required off-diagonal elements 
{ipfAj <p) corresponding to the special choices (A22) of ip{x) and (p{x). The terms involving /2 do not contribute to 
the polo terms because of the Dirichlet boundary condition satisfied by the tost function ip from Eq. (A22). This 
entails that there are no such contributions from the diagonal elements {ip,Ajip), and that there analogs from the 
difference [ij) + (j),Aj[ip + </>]) — {4>,Aj4>) cancel. 



The coefficients /o(l) and /4(1) can be expressed in and 
terms of the single integrals j,j,{m) and jaim) defined by 
Eqs. (47) and (48) as 



/°° df 
dtfo{t)= j^-^fom , (A43) 



/o(l)=<o<^" 



2/2 



and 



/4(1)=<0^'-'"/' 



8-m 



ja{m) 



(A40) 



(A41) 



The two integrals involving /o(i) 

"^dt 



both yield contributions ^ e ^. Let us rewrite them 
in terms of the integrals Ju{m) and ji{m) defined by 
Eqs. (43), (44), and (76). 
To this end, we rewrite Ii as 

I, = lim {Bs + /o(0) In 6), Bs = [ ^ /o(t) . (A44) 



Ii^l f (/o(i)-/o(0)) (A42) 
I 

Using the scaling form (10) of the free bulk propagator and performing the angular integrations, one obtains 
Bs <T™/2 



m >->3-m/2 



-I 



Is Jo Jo ^ ^ 4u;W4(l+^2)(i2 + ,2)2 

Jv(5 



dw 



l/4\ -,00 



2-m/2 



8^1-m/4 (1 + ^2)2-m/4 J^^^^^^ y2-m/i (1 + j.2) _ ^2 ' 



(A45) 
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where 



y{t,r) 



l + r2 



(A46) 



is the new integration variable to which we transformed to get the second Hne of Eq. (A45). To extract the singularity 
of the integral for small 5, we subtract and add '^m,d* (0) in the numerator of the y integral. In the part of the integral 
produced by the term ^rn,d* — ^m,d*(0), we can safely replace the lower integration limit by its 5 — > limit 

2/(0, r). In the double integral J^dr JJ^ ^^dy we reverse the order of integrations whereby it becomes J^dy J^^dr with 
ry = \J y / (1 — y). Upon changing from r to the integration variable x = r/vy, we can perform the x integration. In 
the w integral of the resulting expression, one easily recognizes the scaling function Qmiw^^^) given in Eq.(45) upon 
recalling Eq. (D.l) of Ref. 15. The result of the remaining w integration therefore is proportional to the integral jui'm) 
whose definition is recalled in Eq. (44). 

The part of the integral associated with the constant ^m,d* (0) in the numerator of Eq. (A45) yields terms that 
behave as ~ In^ and approach a finite limit as S 0, respectively. They can be computed in a straightforward 
fashion. Upon gathering all contributions to Bs, substituting them into Eq. (A44), and taking into account that 
/o(0) = (T~™/2 Q, we can perform the required limit S ^ 0. The result is 



./oft) : 

/o(0) . 



J„(m) - 1 - ln2 . 



(A47) 



The integral I2 can be transformed in a similar manner to show that it can be expressed in terms of the integral 
ji (m) introduced in Eq. (76) as 



I ° -m/2 _ / 



foit) . , . 



(A48) 



Combining the above results finally yields 

b b b s 

0--<^^--{® = 0--^^-f® - O-^^-f® + O-^^-f® - O-^^-f® (A49) 



with 



O- 

Zl 



n + 2 ^2/<^^l (PaPa) 



i</>(w) 



24e (8 - m) 



12 (8 -m) 96m (m + 2) 

O(e0)la— -h®, 
J ^1 



(A50) 



O-- 

-21 



O-- 

-21 



(A51) 



a-^-f® 
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-ji(m)+0(e°) 



o~ 

-21 



-h® : 



(A52) 



and 



72 



e (8 - m) 



+ 0(e°) 



O- 

-21 



-h® 



(A53) 



where it should be noticed that the line-symmetry factors of the graphs on the right-hand side of Eq. ( A49) are 1 /2 
(second and third graph) and 1/3! (all others), depending on whether they have two or three equivalent lines. 

We can now substitute the above results into Eq. (73) for G'^'^\p; z) and express the bare variables U and a 
in terms of renormalized ones using the expressions for the bulk renormalization factors Z^f,, Z„, and given in 
Eqs. (40), (41), and (48) of Ref. 15. Upon determining the renormalization factor -^1,00 such that the dimensional 
poles are minimally subtracted, one easily arrives at the series expansion (74). 
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APPENDIX B: RENORMALIZATION OF g£^\p) TO TWO-LOOP ORDER 



In this appendix we consider the renormalization of the two-point function g£'^' (p) for p = r = 0, which requires 
the additive surface counterterm cx C2(?i, e) according to Eq. (73). At two-loop order the sole contribution to this 
counterterm comes from the diagram ®-|— fiiv-]-® . To determine its pole terms one can proceed as in Appendix A, 
choosing the test functions as 



(Bl) 



A straightforward calculation shows that the graph's pole terms in question — i. e., those with no dependence on 
Pjj and a Pa dependence of the form e^^ PaPa — can be written as 



C/2 n + 2 a^/'^PaPa 



8me 9 



f2 



pi POO 

2/2(1) +3y^ dtf2{t) -sj^ dth 



it) 



(B2) 



CT=1 



where U and f2{t) are defined through Eqs. (A9) and (A30), respectively. Upon introducing the coefficient 



4mF2_o r(|-f)r2(i + f) ' 



the two integrals can be written as 

F, " 



p-^ pea 

dtf2{t) = bml^'\m) , -^J^dt /2(t)U=i = bml^^\m) 

7(1) (m) = / dv t;-+i {v) / dy y^+' ^l^^, {vy) f (^) {y^) 
Jo Jo 

pOC pOO 

I^'\rn)=j^ dvv"'+'^m,d'{v)J^ dyy"^+'^l^,,{vy)¥^\y'), 

yt V""-'^/' 



with 



and 



where 



and 



2^/y Vl-y 



(6-m)/4 ^1 ^(2-m)/4 

/ dt , 

Jo \/r^ 



(6-m)/4 



dt 



^(2-m)/4 



VT+I 



2vy V?/-i 

The latter functions are expressible in terms of hypcrgeometric functions as 

^yl^m/A „ y)(m-6)/4 3 ^ 



yt 



(m-6)/4 



2r(2-m/4) 



3 m 3 ^ 
2"T'2^T' ~ T' y- 1 



and 



/7Fr(l-m/4) 



2(y- l)i/2r[(6-m)/4)] 



11m y 



2' 2' 4 1 



+ 



yl-W4(y_l)(m-6)/4 3 ^ 



2V^ 

'3 m 3 TO 



V2 4 ' 2 4 ' 4 ' 



i) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



(B9) 



(BIO) 
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From the above results, Eq. (75) for (72(7/, e) follows at once. 

As a consistency check we have explicitly verified to two-loop order that all poles of gS'^^ (p; ct = t = 0) get canceled 
in the renormalized function G^x>'^ln{p; a = r = 0) for our choice of bulk and surface renormalization functions. For 
interested readers wishing to check this, we here give the Laurent expansions of the corresponding Feynman diagrams 
to the required orders in e and U. They read 



Of-'— --H8 = 



n - 



_2 [/ 
~ 2" 



[l-elnkp + 0{e^)] , 



(Bll) 



-1« 



n + 2 



4 ^ 



■kp) 



1 l-41n(2/tp) 



2e2 



+ 



4e 



+ 0(e°) 



(B12) 



and 



(^) ' T 2i^ { ' + ^ - ' + ^(^')] } 



®-|-f?->H-® = 



6 3 i AkpC 8-m 32m{m + 2)kp e 



V(TPaPa , 
6 Or, 



-:^^-/«(m)+/(2)(m) 
3 6-m ^ ^ ^ ^ 



+ 



3 
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- + Ju{m) - ji{m) - 2 ln(2«;p) 



+ 0(6°) 



(B13) 



(B14) 



where J2,3(w) denotes a particular one of the integrals (49). 



APPENDIX C: ONE-LOOP CALCULATION FOR GENERAL VALUES OF A 



In this appendix we describe the calculation of the renormalization functions Zi{u, A, e), Zc{u, A, e), and P\{u, A, e), 
defined by Eqs. (18), to one-loop order. For this purpose, it is sufficient to consider the cumulants &-^'^\p;zi) and 
dG^'^'°\p;zi)/dc for c = and generic A > 0. Let us indicate that we have set c = by a subscript zero. The free 
propagator (7) in this case reduces to 



Go{p; 2:1,22) 



2kr: 



kp + XpaPa 



Consider first Gq^'^\p; zi). Up to one-loop order, we have 



■ii n + 2 



with 



^ , , f d^-'p k. 
Here ii(A, e; m) is defined by 



kp + XpocPa 2 3 

d'''~^p kp — XpaPa er"^"-^^ 



f 

Jo 



dz Go{p; zi, z) Go(p; z, 0) Go{x, x) + 0{u^) 



+ X PaPa 



ri^-V2 F r(2-e) sin(e7r/2) ^ 
111 Act ' ,e\m)tm,e tzzja z 



^l(ACT-V^e;m).Gl(e)CTW^/^| 



XpaPa e 



+ XpaPcx 



, ii(0, e;m) = 1 . 



(CI) 



(C2) 



(C3) 



(C4) 



where the normalization constant Ci is fixed by the specified value of ii(0, e;m). Thus the A dependence of the 
generalized function Go(a3,a;) is entirely contained in the prefactor ii(A ct^^/^, e; m). 
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After performing the angular integrations in Eq. (C4), we transform from the radial integration variables 
P = {ppppY^'^ and k = {(ipaPaY^^ to k = kp and k, using dkj^ dp g{k,p) = dnj^^dk k {k^ — 
k'^)~^^^ g[k,^K'^ — fe^). The k integration then can be performed, and one easily convinces oneself that ii{\,e;ni) 
can be written in the form (24). 

We can now employ Eq. (A6) to perform the z integration in Eq. (C2). This yields 



12 



(Act 



-1/2., 



1 



1 



2XpaPa 
+ XpaPo 



+ 0(eO) 



Kp + XpaPo 



0{u 



(C5) 



To determine Zi, we can set pa — Q and determine the 0{u) term of Zi{u,\^e) from the condition that the pole 
of G^q'^\p\ zi) cancels in G^^'^J^{p\zi). This yields the result (21). For pa 0, the pole imphed by the first term 
of Eq. (C5) in parentheses cancels upon multiplication oi ' . Requiring that the remaining one implied by the 
second term in parentheses gets absorbed by the renormalization of A according to Eq. (18) yields Eq. (23). 

Turning to the function [^(^(^'"^(p; zi, Z2)/9c]o, we note that the derivative —d/dc generates an insertion of the 
surface operator J ^dA<p'^ /2, which we depict as . Hence we have 



dc 



?G(2'0)(p;.r,.2) 



c=0 



Zl 



Z2 



a 



Zl 



-O + {zi ^ Z2) 



Zl 



Zl 



Z2 



(C6) 



The Feynman integral of the graph inside the brackets is nothing but G(p; zi, 0) times the term ex u of Eq. (C2), with 
the right external point taken at Z2- The remaining last graph's loop involves the integral 



' {kp + Xpc^Poc) 
where 12 is the counterpart of ii defined by 

Z2(Aa-i/2^e;m) =C2(e)(7"/4 



Frr. 



4r(l-cj sm(r^/2) . • 



m/i 



(2^) (Kp + ApaP„)' 



, 12(0, e;TO) = 1 . 



(C7) 



(C8) 



Proceeding as in the case of ii, one arrives at the form (25). 

A combination of the above results yields the Laurent expansion 



^G(2'0)(p;zi,Z2) 



1 



Un + 2 



c=0 



»i(Aa-V2) 1 



I e 6 

-2i2{\a-^'^;m)+0{e'') 



_ 2\paPa 
kp + XpaPa 

g-Kp (21+22) 



(C9) 



{kp + XpaPa) 

from which Zc can be determined in a straightforward fashion by requiring that the poles of the func- 
tion (C9) are minimally absorbed through Zc so that the renormalized quantity [9Gren'^(p; zi, Z2)/9c]c=o = 
Zc Z^^ [9G^^'°)(p; Zi,Z2)/dc]o = [dGren\p; Zi,Z2)/dc]o becomes uv finite. The result is given in Eq. (22). 

Using Mathematica,^^ the integrals ii(A;m) = ii(A, e = 0;m) and i2(A;m) = i2(A, e = 0;m) introduced in 
Eqs. (24) and (25) can be computed for < m < 6 and expressed in terms of hypergeometric functions. One obtains 

2tt [1 - (1 - A2)(2-™)/41 



ii(A;m) 



m 



( m 5 ,r,\ (mi 
2Fr(l,l + -;-;A^)+2Fr(l,-;-; 



m 3 



A^ 



and 



(Vft 1 



A cos(m7r/4) i3[m/4, (6 - m)/4] 
^A(l-A2)-(^+'")/^ 



cos(m7r/4) i?[(2 - m)/4, m/4] 



(CIO) 



(Cll) 
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For m — > 0, m = 2, m = 4, and m 6, these expressions simplify to 



n(0+,m) = 1, (C12) 

ii(A;2) = -1 + ^ ln(l + A) , (C13) 
A 

TT - A 2arccosA ,^ 

zi(A;4) = — ^ , (C14) 



-A 
l + A ' 



n(A;6-) = (C15) 



and 



i2(0+;m) = 1, (C16) 

i2(A;2) = (C17) 

. /, ,^ 1 — AarccosA 

i2(A;4) = , (C18) 



1 - A2 
1 



*2(A;6-) = 7:^^, (C19) 



respectively. 



APPENDIX D: CALCULATION OF "^'(0) 

In this appendix we wish to determine the renormalization (18) of A to order for A = and compute the residuum 
'"^^(0) introduced by Eq. (20). 

To this end, we consider the function Gq^'^\p; zi) = G^^2^{p',zi) for c = c = and A = 0. As claimed, the bare 
variable A does not vanish at order if A = 0. Hence we must keep a nonzero A in our calculation. Up to two-loop 
order, the Feynman graph expansion of G'q'^\p;zi) reads 

G^^'^\p;zi) = o ®+o Q ® + O 00 ^ + p__8 8 

Zi Z\ Z\ Z\ 

+ O — ® + O(m^) , (Dl) 

z\ 

where O O here represents the free propagator (7) for c = f = 0. As before, the crossed circle denotes a point 

z z 

on the surface. Thus the first graph simply becomes 

-® = » • (U^J 



Z\ Kp + XpaPa 

In Appendix C we saw that at one-loop order the renormalization of A remains multiplicative. Unless A ^ 0, the first 
two two-loop graphs in Eq. (Dl) obviously cannot produce primitive momentum dependent uv singularities of a form 
corresponding to the renormalization of the boundary operator X^aC^"?^)^ because the (tadpole) graphs one obtains 
through amputation of the external lines arc independent of the momentum p. This momentum independence holds, 
of course, also if A ^ 0. In this case, uv singularities requiring a multiplicative renormalization of A are produced by 
the pole ~ S'{z) of the amputated one-loop tadpole graph. As we explicitly showed in Appendix C, this follows quite 
simply by computing the action of this distribution on the external legs. This does not happen for A = as long 
as A may be set to zero in the free propagator so that it obeys Neumann boundary conditions. The upshot of these 
considerations is that the first graph contributing to the (non-multiplicative) r enor malization of A is the right-most 
two- loop graph in Eq. (Dl), involving the momentum dependent subgraph ■ 

Since we know that A = O(u^) if A = 0, we can set A = in all but the first graph of Eq. (Dl). Thus the lines 
of all those graphs correspond to the free Neumann propagator, which in analogy to Eq. (Al) can be written as 
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Gn = Gb + Gs- Using this decomposition and the graphical conventions introduced in Eqs. (A2) and (A49), we may 
split the two-loop graph in question as 

s 

+ O (J^ ® , (D3) 

Each term on the right-hand side can be computed in the manner described in Appendix D. One chooses the functions 
and 'ip as^® 

utilizes Eqs. (A23)-(A25), (A29), and (A37)-(A39), and finally computes the Fourier transform / d'^-Vi2 • • . e'^'^'i^ 
Only those parts of the quantities Ij that are proportional to / d'^~^r (fi{r, 0) dadaf{r, 0) contribute to the renor- 
malization of A at order if A = 0. They are given by 

lo = jd''-'rip{r,0)dadMr,0) + ... , 

/i = dtf2{t) Jd''-^r^{r,0)dc,dMr,0) + .-- , 

h = -^J^dth{t) Jd''-'r^{r,0)d„dMr,0) + -- - , 

h = ^ jd''-'r^{r,0)d^dMr,0) + ... , (D5) 

where the ellipses stand for other types of terms. Upon going over to the pz representation, we obtain the result 

. . r°r2 „ I o n-™/2 r°° rfi 



The explicitly displayed pole term must be absorbed by the rcnormalization of A according to Eq. (18). 
The expansion of the zero-loop graph o ® to first order in A produces a contribution of the form 

P\{u, A=0, e) (7^/^ [d&^'^'> Its u^/e part must cancel the above pole. This is the case if we make the choice 

.-)(0) = 1^ rdtm = 6„ [7(^)(m) + /(2)(m)] , (D7) 

where I^^^ (m) and 7^^^ (m) are the functions introduced in Appendix B [cf. Eqs. (B4)-(B6)] . 

APPENDIX E: NUMERICAL EVALUATION OF INTEGRALS 

In this appendix we briefly explain how the numerical values of the integrals J2,3(m), ji{m), l[^\m), and I^'^\m) 
presented in Table I were obtained. 

1. Numerical calculation of the integrals J2,3(m) and ji{m) 

From its definition (49) it is clear that J2,3(?Ti) is a single integral of a similar form dvf{v) as the quantities 
j,f,{m), jcr{m), jp{m), and j„(m), which were previously introduced and numerically evaluated in Ref. 15. The same 
applies to the integral ji{m) if we utilize its representation (77) in terms of flm,d'- Unless m takes the special values 
2 or 6, the integrands f{v) involve differences of hypergeometric functions that grow exponentially as v ^ 00, while 
the f{v) themselves decay as inverse powers of v (cf. Appendix E of Ref. 15). 

To cope with these difficulties, we proceeded as follows. Summing up the Taylor series expansions of the scaling 
functions ^m,d' (u) and ^m.d* (v) appearing in the integrands f{v) gave reliable results as long as v was not too large. 
Typically, this worked for all v < vq up to vq ~ 10.7. For larger values of m, even larger choices of vq were possible. 
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We therefore worked with vq = 10 for both functions d* (w) and f^m d. (u), and all m. Depending on whether v 
was smaller or larger than this value of vq, we relied on the Taylor summation and the asymptotic expansions of the 
functions $m,d*(u) and ^m,d'{v) (see below). 

In order to gain precision and speed, it proved useful to compute the required Taylor-series expansion coefficients 
in a recursive fashion. Starting from Eq. (12), one can easily show that the function ^m,d' can be written as 



— ^ E^4-V64) -^r^.v- 



(El) 



where Sfe(a;) and s{x) satisfy the recursion relations 



8xsfe(x) . 4a;sfc(a;) 

^''+^(^^= (l + 2fc)(2 + m + 4fc) ' ^'-+^(^^= (l + fe)(4 + m + 4fc) ' -o(a:) = so(x) = 1 . (E2) 



Prom Eq. (79) one easily derives the following analogous representation of 



"™'''*^^^-8r(l + m/4)^ 2-m-2 + 4fc-r(i + ^)'' ^ m - 4 + 4;fc + ^^^^ ^^^^ 



with 



and 



^-(^) (TTwSt^ ' ^-(^) = (rTWT^' W = ^o(.) = i, (E4) 

8(m-2)r(l + f) 32 (m - 2) sin (m7r/4) (TO-4)r(| + f)' ^ ^ 

The apparent singularities at m = 2 and m = 4 actually cancel, so that 

. ^ / k + HCE~l) + k Hv/2) for m = 2 . 

"^""^ I ^7ri/2i;2 + 7r-i/2[l-CB+ln(2t;-2)] for m = 4 . ^ ^ 

In order to determine the values of the scaling functions for v > vq, we utilized their asymptotic expansions. Those 
of the functions ^m,d' can be found in Eqs. (A5) and (A6) of Ref. 15; those of the fim^d* (v) can be shown to read 

O M- ^\,rn-SY-f 64^ r(§ + fc) 

n^Mv)--^v ^^(^-^j (f-fc-2)r(f-i-fc)- 

Note that if m = 2 and m = 6, this asymptotic series truncates at zeroth and first order, respectively, just as its 
analog for ^m.d* {'>^) does for these choices of m. 

Keeping the first three terms of the asymptotic expansions turned out to be sufficient. At v = vq = 10, the Taylor 
series of both ^m,d* {v) and $m,d* {v) agreed with their respective (truncated) asymptotic series to within 10~^ to 
10"'' percent. 

Employing the above methods to compute the values of the integrands f{v), the required integrals could be de- 
termined by straightforward numerical integration, giving the results for J2,3(to) and ji{m) presented in Table I. 
The numerical values found for m = 2 and m = 6 in this manner agree with the analytical results to the number 
of decimal digits retained or better. We also computed ji{m.) by numerical evaluation of the double integral (76), 
obtaining results in conformity with those obtained by the above method. 



2. Calculation of l'-^\m), I^^\m), and related 
quantities 

(2 —1) 

The renormalization hmctions C2(w,, e) and ' (0) 
involve the double integrals l'^^\m) and I^'^\m) defined 



by Eqs. (B5) and (B6). In the special cases to = 2 and 
TO 6, the required integrations can again be done an- 
alytically. The functions 4-^^ and [cf. Eqs. (B9) 
and (BIO)] appearing in the integrands of I^^\m) and 
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/(^^ (m) reduce to the simple expressions 



(2), 



1-y 



1/21 



(E8) 



and 



^i^Ky) = = g^T;^ [y-'^' + o(6 - m)] , (eq) 

respectively. Upon substituting these results together 
with the simplified expressions for the scaling functions 
*1^2,5(^') and $6, 7(1') into Eqs. (B5) and (B6), the integra- 
tions can be performed with the aid of Mathematica®-^ 
to obtain the results (84)-(87). From them Eqs. (59) and 

(60) for Pi^'^^\0) follow by insertion into Eq. (D7). 

For other values of m, we performed the dou- 
ble integrals (B5) and (B6) numerically, employing 
Mathematica.®^ In contrast to ^m,d'{y), the functions 

^rn (y) and '^ir} (y) are not differences of functions that 
grow exponentially as y ^ 00. Therefore, their calcula- 
tion poses no problem for MATHEMATICA even for largo 
y. In the case of iv) (and its powers appearing in 
the integrands), we simply let MATHEMATICA evaluate 
its representation (12) in terms of hypergeometric func- 
tions if 2/ < yo ~ 8.5, and utilized its asymptotic large-?/ 
expansion for y > yo. In each case we made sure that 
the values of the integrands the two methods yielded at 
the matching point y = yo were sufficiently close so that 
their difference could safely be neglected, given the lim- 
ited number of decimal digits of our final numerical re- 
sults. 

The integrands of I^^\m) and I^'^\m) have integrable 
singularities at 2/ = 1 and vary considerably in some 
parts of the two-dimensional integration regime. How- 
ever, Mathematica's algorithm for numerical integra- 
tion requires that the integrand does not vary too much. 
Letting us guide by two-dimensional plots of the inte- 
grands, wo therefore divided the integration regimes into 
suitable rectangular subregions such that the algorithm 
was able to perform the integration over each of them 
without running into precision problems. We then added 



up the contributions from these rectangles to determine 
the integrals l'-^-' (m) and J*^^^ (m). The resulting numbers 
can be found in Table I. 

The numerical values Mathematica produces for 
l'"'^-' (m) and I'^^-' (m) in the special cases m = 2 and m = 6 
agree with the analytical results (84)~(87) to within 10~® 
percent. However, since Mathematica automatically 
makes use of the fact that scaling functions such as ^m,d* 
simplify considerably if m = 2 and m = 6, we must cau- 
tion not to conclude that the same high numerical preci- 
sion applies to our numerical results for other choices of 
m. 

Let us also note that in the case n 
^m iy) simplify to 



4, *2^(2/) and 



and 



*f (2/) 



^t>{y) = {l-y)-''^K[yl{y-l)] 



(ElO) 



(y_ 1)1/2 _yl/2 



(y- 1)1/2 +2/1/2 



(y_ 1)1/2 +yl/2 

(Ell) 

where K[y) denotes a complete elliptic integral of the 
first kind. 

A final remark concerns the sum /(i^(m)-|-/(^^(m) ap- 
f 2 — 11 

pearmg m PI '{0). From Eq. (B4) we see that it is 
proportional to J^dtf2{t). The integration over t is 
straightforward, producing a free propagator in rf* — 1 
dimensions. The result translates into 

/(i)(m)+/(2)(m) 

4 Jo Jq 

x{R' + l)-'^l^,,[g{R' + l)-'/^]. (E12) 

Since the integrand on the right-hand side has less struc- 
ture than those of the individual integrals I^^^m) and 

/•^^^(m), the form (E12) of the sum lends itself more eas- 
ily to numerical evaluation via Mathematica. Our re- 
f 2 — 11 

suits for PJ; ' ''(0) presented in Eq. (61) were obtained 
in this fashion. 
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